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Abstract 



Qh We study the transport properties of topological insulators, encoding them in a 

generating functional of gauge and gravitational sources. Much of our focus is on 
the simple example of a free massive Dirac fermion, the so-called Chern insulator, 
especially in 2+1 dimensions. In such cases, when parity and time-reversal symmetry 
{SI are broken, it is necessary to consider the gravitational sources to include a frame 

and an independent spin connection with torsion. In 2+1 dimensions, the simplest 
parity-odd response is the Hall viscosity. We compute the Hall viscosity of the Chern 
insulator using a careful regularization scheme, and find that although the Hall viscosity 
fNj is generally divergent, the difference in Hall viscosities of distinct topological phases 

is well-defined and determined by the mass gap. Furthermore, on a 1+1-dimensional 
edge between topological phases, the jump in the Hall viscosity across the interface is 
encoded, through familiar anomaly inflow mechanisms, in the structure of anomalies. In 
particular, we find new torsional contributions to the covariant diffeomorphism anomaly 
in 1+1 dimensions. Including parity-even contributions, we find that the renormalized 
generating functionals of the two topological phases differ by a chiral gravity action with 
a negative cosmological constant. This (non-dynamical) chiral gravity action and the 
corresponding physics of the interface theory is reminiscent of well-known properties 
of dynamical holographic gravitational systems. Finally, we consider some properties 
of spectral flow of the edge theory driven by torsional dislocations. 
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1 Introduction 



Quantum field theory anomalies imply that symmetries that were present in the classical 
Lagrangian are broken due to quantum effects. While at one time they might have been 
thought of as a sickness of certain field theories, anomalies lie at the heart of some of 
the most fundamental physical phenomena in real materials. The canonical example is 
the integer quantum Hall effect (IQHE) where a 2+1-dimensional electron gas in a large, 
uniform magnetic field exhibits a Hall conductance which is quantized in units of e 2 /h when 
the chemical potential lies in a Landau level gap (and has been measured to be quantized 
up to 10 significant digits). The precise quantization arises from the connection between 
the Hall conductance and a topological invariant of 2+1-d electron systems called the first 
Chern number C\. Since C\ is a topological quantity which is determined by the ground 
state, it is not affected when the system is perturbed continuously, and is insensitive to the 
microscopic details of the sample as long as the bulk energy-gap is not destroyed. Thus, 
response coefficients that are determined by topological invariants are the most universal 
features of gapped systems. 

For all understood topological response coefficients there is a complementary way to view 
the quantization by studying the properties of the gapless, fermionic modes that lie on the 
boundary of the system. There is a deep connection between topological transport in the bulk 
of a gapped material (say in 2+1-d) and field theory anomalies that are present for the (say 
1+1-d) gapless boundary states[l, 2]. The connection between anomalous currents, topology, 
and index theorems underlies some of the most beautiful transport phenomena that have been 
predicted, and in some cases observed in real materials. For the IQHE this bulk-boundary 
correspondence connects the bulk Hall transport to the spectral flow of the boundary chiral 
modes due to the chiral anomaly. The edge anomaly provides a complementary picture of 
the origin of the Hall conductance quantization which is commonly known as Laughlin's 
gauge argument (though it was not originally written in terms of anomalies) [3]. 

While most anomalies connected with charge and spin currents are well understood, the 
anomalous thermal, and visco-elastic responses (VE) are not. The thermal and VE re- 
sponses lie at the intersection between geometry, topology, and quantum field theory as they 
are usually represented as topological phenomena associated to geometric deformations of a 
field theory. One example of such a novel effect is a dissipationless, electronic viscosity re- 
sponse in the 2+1-d topological Chern insulator with broken time-reversal symmetry[4, 5, 6]. 
While the ordinary shear viscosity generates a frictional force tangent to fluid motion, the 
dissipationless viscosity produces a perpendicular force (see Figure 1)[4, 7]. This viscosity is 
not clearly understood except in some special cases including the integer and fractional QHE 
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Figure 1: Fluid mechanics illustration of the viscous forces. A counter-clockwise rotating 
solid cylinder immersed in 2d liquid droplet with (a) non-zero shear viscosity (b) non-zero 
dissipationless viscosity. Note that the resulting forces (arrows outside cylinder) are tangent 
and perpendicular to the cylinder motion (arrows inside cylinder) respectively. The shear 
viscosity impedes the cylinder while the dissipationless viscosity pushes fluid toward or away 
from the cylinder depending on the rotation direction. 

with rotation[8, 9] and translation invariance[10], and chiral superconductors [1 ']. However, 
all of these models share the feature that they are Galilean invariant, and in relativistic 
systems, or lattice models with broken continuous translation symmetry, it is not clear if 
the topological viscosity is quantized, or even well-defined (for the lattice case) [6]. This is 
unusual as one would expect that it should be quantized like all of the other examples of 
topological response coefficients, such as the quantized Hall conductance (which is simulta- 
neously present in the 2+1-d Chern insulator phase) [12]. 

In this article we will address these issues by constructing an explicit bulk-boundary corre- 
spondence which allows us to understand the anomaly mechanism associated to the topo- 
logical viscosity. The interplay of the topological response with the geometric deformations 
of the system makes this problem more subtle than previous known examples of topological 
responses, because, while topology does not care about the details of a shape, geometry does. 
The bulk-boundary correspondence for the viscosity response is completely unknown and, 
as we will indicate below, must have a different physical origin than the, say, chiral anomaly. 
The model we will focus on for most of this work is the massive Dirac model. This model 
represents the low-energy physics of topological insulators in various dimensions, and with 
various symmetries [] ;]. This model responds quite differently to geometric perturbations 
than typical non-relativistic electrons (i.e. systems with small spin-orbit coupling). To il- 
lustrate the underlying premise, we first note that conventional non-relativistic electrons in 
a crystal are described by the Schrodinger equation at low-energy, and are only elastically 
influenced by the stretching of bonds that is captured by the strain tensor[l []. However, 
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spin-orbit coupled electrons described, for example, by the Dirac equation at low-energy, 
are also aware of the local orbital orientation, which is not contained in the strain tensor. 
Instead the Dirac model couples to geometric perturbations via a local "frame field" that 
we will introduce below. This additional sensitivity generates physical responses to shear- 
ing, twisting, and compressing/stretching that are not found in weakly spin-orbit coupled 
systems. These phenomena are the focus of our work and are connected with the idea of 
geometric torsion as we will discuss. 

Our paper is organized as follows: In section 2, we introduce some basic concepts of geome- 
try and elasticity, with a special emphasis on torsion, from a condensed matter perspective 
that are relevant to our later discussions. We will follow this up with a more mathemati- 
cally precise description in section 3, from the point of view of general relativity and high 
energy physics. In section 4, we introduce some basic aspects of fermions in the presence of 
background gauge and gravitational fields, again focussing on the role of torsion. Through 
sections 3 and 4, we will also set up notation that will be used in the rest of the paper. 
In section 5, we calculate and carefully regularize the Hall viscosity for the Dirac model in 
2+1 dimensions, working about a flat background. In section 6, we will then proceed to 
compute the full effective action in the large mass limit on a generic background, and show 
the emergence of the chiral gravity action in the non-trivial topological insulator phase. One 
of the main ideas in sections 5 and 6 will be, that the differences in transport coefficients 
between different phases are physically meaningful. We will explore this further in sections 7 
and 8, where we study the parity-odd transport properties in the context of anomalies due to 
chiral edge states localized on the interface separating a non-trivial phase from a trivial one. 
Finally, we will study Hall viscosity from the point of view of the interface Hamiltonian spec- 
tral flow in section 9, and discuss possible mechanisms for the corresponding bulk-boundary 
momentum transfer. 



2 Informal Preliminaries 



Before we move on to a more precise description with which high-energy theorists will be 
more comfortable, we try to informally introduce the necessary background material for a 
condensed-matter audience using the language of elasticity theory. Conventional elasticity 
theory is one of the foundational underpinnings of solid state physics as it contains within it 
the physics of the lattice structure, including, for example, phonon fluctuations away from 
the ordered reference state. At a given time, one characterizes an elastic medium via a 
displacement field u(x n ) which gives the vector displacement of a lattice site n, away from 
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Figure 2: (a) Reference state (hollow circles) and displaced state (solid circles) for an elastic 
medium. Displacement vectors for each site n are denoted by u(x n ). Zoom-in shows frame 
field vectors e 1; e 2 in the reference state (aligned to crystal x, y-axes) and the displaced state 
(rotated with respect to crystal axes), (b) Edge dislocation representing the fundamental 
torsion lattice defect. An electron traveling the thick line surrounding the dislocation will 
be translated with respect to the same path in the reference state that does not enclose a 
dislocation. The Burgers vector is in the y-direction. (c) Disclination represented by a single 
triangular plaquette in a square lattice crystal. Gives rise to curvature i.e. objects that 
travel around a disclination are rotated with respect to the reference-state path. 



the position x n of a given reference state (note that we will take the continuum limit where 
n becomes a continuous label and thus x n becomes a continuous coordinate yielding a field 
u(x)) . If every lattice point is displaced by the same amount then the crystal has just been 
globally translated and does not feel any internal stress. However, if the displacements of 
lattice sites are not identical, the material will respond by generating a stress (momentum- 
current density) 

T ij = A ijH u u + V ljki u ke , u M = l/2(d k u e + d e u k ) (1) 

where repeated indices are always summed, T u is the stress tensor (momentum current den- 
sity), A y is the elasticity tensor which relates stress to the strain u k i {i.e. a generalization 
of Hooke's law) and rfi is the viscosity tensor relating stress to the strain rate/velocity 
gradient itke (i.e. a velocity dependent frictional force). See Figure 2a for an illustration of 
a lattice elastic medium and a displacement field. 

A non-zero strain tensor indicates that the (spatial) geometry of the elastic medium has 
been distorted. The geometric characterization of the lattice is contained in the metric 
tensor which determines the distance between lattice points. In the ordered reference state 
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shown in Figure 2a the metric tensor is just = 5y which implies that distances between 
sites are calculated in the usual Euclidean way. When the material is strained, the spatial 
metric tensor is modified to become gij = 8{j + 2iijj[14], which is what is meant when we say 
the geometry is deformed. Static lattice deformations affect the electronic behavior since the 
bonds are deformed. For electrons described by the Schrodinger equation at low-energy, the 
Hamiltonian is modified to become (to linear order in strain) 

2m 2m 2m v 1 2m v n "m y ' 

where g l ^(x) is the inverse of the metric tensor which depends on position via the contribution 
of the strain tensor. Thus, depending on the spatial profile of the strain, the electron 
spectrum can be drastically modified. 

While the strain/metric based elasticity theory is quite successful, it is not general enough 
to model all of the electronic structure effects arising from the coupling of materials with 
spin-orbit coupling to geometric deformations. What is needed is a more fundamental field: 
the frame field e a in d spatial dimensions where a = 1,2, ... d labels each vector of the frame 
(with components e^). The frame- field is a set of d vectors residing on each lattice site, 
and heuristically encodes the local bond stretching (through the vector lengths) and the 
local orbital orientation (through their relative angles on each site). As we will see later, 
in many instances it is more natural to consider the co-frame field e a which is a local basis 
of 1-forms that are dual to the vectors e b (i.e. they satisfy e a (e_ b ) = 5%). For the reference 
state shown in Figure 2a the reference frame fields are orthonormal vectors which are aligned 
with the crystal axes. The distances between lattice sites, i.e. the (inverse) metric tensor is 
determined from the frame fields via g^(x) = S e* (x)e^(a;)[15]. It is easy to see that if the 
frame fields are orthonormal at each site then g % i = 6^ as expected. The key relationship 
between the metric and the frame is that we can locally rotate the frame at each site by any 
SO(d) rotation matrix R and we get the same metric back: 

& = 5 ab (R c a (xU(x))(R d b (x)^(x)) = R c a (x)R d a (x)e^ d = 5 cd £e d = g* (3) 

since RR T = I. This implies that an elasticity theory determined completely from the met- 
ric does not capture local orbital deformations since each different local orbital orientation 
yields the same metric tensor. However, electrons with SOC propagating in a lattice will be 
sensitive to the local orbital orientation, which is exactly why a frame field must be intro- 
duced to couple these materials to geometric perturbations. This modification to elasticity 
theory is closely related to so-called micro-polar or 'Cosserat' elasticity[16, 17]. 

At this point it is useful to explicitly show how the frame field enters spin-orbit coupled 
Hamiltonians. The low-energy description of two such systems are given by the Dirac 
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Hamiltonian (which represents, for example, topological insulators) [18, 19] and the Lut- 
tinger Hamiltonian (which represents, for example, the upper-most valence bands of III-V 
semiconductors) [20, 21]: 

Hd = w^^elr + mr (4) 

La 

Hl = ^^ + a(p^)fr^ (5) 

for Dirac matrices T a , spin-3/2 matrices S a , and parameters v, m, a. Hence, the prescription 
is to replace terms of the form piM 1 for a matrix M l , which arise naturally in materials with 
SOC, with J2 a Pie l a M a . Note that for H L , since S a S b ^ 5 ab , the quadratically dispersing 
Luttinger model is indeed affected by the local orbital orientation since it couples to more 
than just the metric tensor. The effects of the frame field are thus not limited to the 
linearly dispersing Dirac equation and affect any coupling between the direction of electron 
propagation pi and the spin/orbital degrees of freedom represented by M l . 

There are two complimentary interpretations of the (co-)frame-field which we will use. The 
first interpretation is in terms of familiar elasticity quantities, namely to first order in the 
displacement field, the co-frame and frame can be expanded as 

1 l + dx^ ~ a a dx t {) 

where diU a = wf is the distortion tensor which is familiar from elasticity theory [14]. The 
quantity wf is effectively the unsymmetrized strain tensor and contains information about 
local rotations through the anti-symmetric combination My = 5i a Wj — 8j a w\. The distortion 
tensor also contains information about dislocations through the line-integral 



wtdx 1 = f du a = -b a (7) 
c Jc 

where b a are the components of the total Burgers vector of the dislocation(s) enclosed within 
the curve C (see Figure 2b for an example) [14]. 

For point-like dislocations in 2d we can write de a = —b a 5^{x) from Stokes' theorem where 
de a is the exterior derivative of the 1-form e a . This formula suggests a second description 
of the e a as a set of d vector potentials. As a comparison, we know that for electrons in 
an electromagnetic vector potential we use the minimal coupling replacement pi — > p\ + qAi 
which shifts the momentum in the Hamiltonian, and we have already mentioned that the 
proper replacement for the frame field is to scale momentum 

Pi -»> Pi£ = PiS* a ~ Pi W i =Pa~ ViWi- (8) 
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Momentum Transfer 




Figure 3: Laughlin gauge argument for torsion: Thought experiment with an insertion of 
torsion flux i.e. a dislocation into cylindrical hole, equivalent to twisting cylinder in the 
y-direction as a function of time. Non-zero dissipationless viscosity causes transfer of p y - 
momentum in the x-direction, i.e. a momentum current perpendicular to time-dependent 
strain. 

Comparing to the electromagnetic case, this shows that each frame-vector yields a vector 
potential that minimally couples to electrons via momentum i. e. the momentum components 
are the charges of these gauge fields. With this interpretation, dislocations are just the 
magnetic fluxes of these vector potentials, and the translation effect of a dislocation is just 
the Aharonov-Bohm effect for the co-frame vector potentials. In general we can construct 
the torsion tensor, which, in the absence of curvature can be chosen to take the simple form 
of a field strength tensor of the co-frame vector potentials Tij a = die" — djef. This has an 
extra index a compared to the electromagnetic version Fij, which labels the particular vector 
potential/co-frame potential. This is how "torsion" naturally enters the discussion, and as 
we can see, it is intimately connected to dislocation density. We also mention that there 
exist other elastic defects such as disclination defects which represent sources of geometric 
curvature and also orbital-twisting defects that can be produced in a strain-free lattice with 
a trivial metric but non-trivial frame (e.g. a torsional skyrmion[22]), examples of the former 
are shown in Fig. 2c. 

With the background theory now set up, we will move on to discuss the current state of the 
field of topological VE response, and some of the open questions which we are proposing to 
study. The first calculation of a topological VE response was the work of Avron et al. which 
showed that a dissipationless viscosity is present in integer quantum Hall states [4, 23, 7]. 
The work was not followed up on until over a decade later when Read showed that, not 
surprisingly, fractional quantum Hall states could also exhibit such a viscosity, and that 
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the response would be quantized if rotation symmetry were preserved[8, 9]. Soon after, 
Haldane showed that rotation symmetry is not a necessary ingredient and that the viscosity 
is related to a fundamental property of an unreconstructed quantum Hall edge: the edge 
dipole moment [10]. For these systems the viscosity, denoted (h, is a quantized multiple of 
h/£ 2 B where £b is the magnetic length. This quantity has units of angular momentum density, 
or momentum per unit length, or dynamic viscosity (force/velocity), and interestingly, it 
depends on a non-universal length scale which varies when the magnetic field is tuned. 
In fact, one even can remain on the same Hall plateau with fixed Hall conductance, and 
tune the field so that the viscosity changes. When rotation symmetry is present, conserved 
angular momentum can be transferred between edges via an applied torque (e.g. due to the 
electric field generated from perpendicular applied flux). The amount of transferred angular 
momentum does not depend on £b, and is given by the quantized multiple of h appearing 
in (h. The same is true of the edge dipole moment, which is also independent of is for 
unreconstructed edges, and is the same universal number multiplying h. 

This quantization emerges quite naturally in the Landau level problem where the quantum 
Hall effect is generated by an external magnetic field. However, the situation is much more 
subtle and complicated when the quantum Hall effect is generated by a topological band 
structure which can naturally furnish multiple length scales. We will focus on this type of 
system to study the impact that a combination of geometry and topology will have in band 
theory. The topological viscosity has been calculated in a (properly regularized) continuum 
model for the Chern insulator, i.e. the massive Dirac Hamiltonian in 2+ld [6]; the details 
of this calculation will be presented in section 5. The regularized value was found to be 
(h = where £ = hv/2m is the length scale induced by the Dirac mass m (with units of 
energy) for a material with a Fermi- velocity (speed of light) v. In relation to the discussion of 
elasticity theory above, the non-zero viscosity coefficient produces a Chern-Simons response 
for the co-frame fields: 

SeffK\ = y J d 2 xdt e^e;d v e b pVab (9) 

where a,b = 0,1,2, and r\ ah = diag[— 1, 1, 1] is the flat-space Minkowski metric. This is 
essentially multiple copies of the conventional Abelian Chern-Simons term, one for each of 
the co- frame fields (including the co-frame in the time direction). As shown in Ref. [6], 
if we calculate the electronic contribution to the stress response T£ = de t( e ) , one finds 
that electron momentum- density is bound at dislocation defects and momentum- current is 
generated perpendicular to any velocity-gradients/strain-rates (see Figure 3a for a picture of 
the latter). This is completely analogous to the charge density bound to magnetic flux and 
charge current produced by electric fields (or time- dependent fluxes) in the quantum Hall 
Chern-Simons response. Additional discussions of the Hall viscosity in topological insulators 



11 



and the relation to lattice deformations and the electron-phonon coupling[ ], a spin-Hall 
viscosity [25], and a Streda-like formula for the viscosity[26] have been carried out. 

The principal issue is that the viscosity does not appear to be quantized, or even universal, 
which is very strange in light of all the previous results on topological responses in topological 
insulators[27, 12]. In this paper we approach these issues by considering the bulk-boundary 
correspondence between momentum-transport in the bulk and anomalous currents in the 
chiral edge states. We develop a picture analogous to the charge response for the quantum 
Hall effect in terms of bulk Chern-Simons response and a boundary chiral-anomaly. To 
address these issues we must develop more precise language than we have used in this more 
informal section; we do this now. In particular, we will introduce a connection to account 
for parallel transport, which will allow us to restore Lorentz covariance and treat torsion 
and curvature on the same footing. As a note, a condensed-matter minded reader might 
first read sections 5, 6, and 9 to get some picture of the physics before tackling the more 
technical, but essential, discussions of the other sections. 



3 Geometry with Torsion 

Gravity is usually described as a theory of metrics, corresponding to a measure of invariant 
distance 

ds 2 = g^dx^dx" (10) 

where x M are local coordinates on a manifold. We can package the information contained in 
the metric (and more in fact) into the components of a co-frame, a local basis of 1-forms 
e a = e a ^dx^ on the manifold. Equivalently, we can regard e a as a local section of the oriented 
co-tangent bundle of the manifold. The metric is related to the components of the 1-forms 
via 

efelVab = 9n» (11) 

where r\ a \, are the components of the Lorentz-invariant Minkowski metric. We will denote 
the dual set of frame vector fields as e a , with e a (e b ) = 5%. To translate (co-)tangent bundle 
data from point to point on the manifold, we need a connection or covariant derivative V. 
Conventionally we write the translation of the frame along a vector field Jf as 

Vxe a = -uj a b {X)e b (12) 

where we have introduced the components of the spin connection u a b, which we regard as 
a set of 1-forms. In a basis of local coordinates this equation can be written as X M V M e" = 
— X^uj^ a b e b u . The spin-connection can be thought of as a non-Abelian gauge field that couples 
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to the rotation and Lorentz transformation generators. Throughout our work, we will make 
one assumption about this connection, which is that it is metric compatible. In metric terms, 
this means that the metric is covariantly constant Vx9 = 0, but using the relationship 
between the metric and the co-frame and the definition (12), it also corresponds to the spin 
connection being valued in the orthogonal group, 1 i.e., u a b = —Ub a (where u a b = r] ac uj c b). 
Under a local change of basis (i.e., a local Lorentz transformation) e a i— > A a be b , the connection 
transforms as 2 

u\ ^ (AuA- 1 -dAA- 1 )\ (13) 

Thus u a b is the 'gauge field' for local Lorentz transformations. The curvature 2-form, or 
field strength, of the connection 

R a b = duj a b + co a c Auj C b. (14) 

transforms linearly 

R\ i y (ARA- 1 )^. (15) 

The components of the curvature 2-form give the Riemann tensor, R cd = \R a b )C d^ a A e b . If 
we denote the covariant derivative acting on (local) Lorentz tensors by D, the torsion 2-form 
is defined as 

T a = De a = de a + u\ A e b (16) 
Torsion also transforms linearly under local Lorentz transformations. 3 

T a h->- A a b T b . (19) 
We write the components of the torsion 2-form as T c = \T c ab e a A e b . 

A very basic property of the connection, is that it satisfies the following translation algebra 

[V a ,V 6 ] = -T c ab V c + R cd . ab J cd (20) 

1 We will work in d = D + 1 spacetime dimensions, so the relevant orthogonal group is SO(l,D). 

2 Note that we are reserving the term 'Lorentz transformation' for these local changes of basis for the 
orthonormal co-frame. These should not be confused with (linear) diffcomorphisms, which are local changes 
of the coordinates. 

3 The Bianchi identities are 

DR a b = dR a b - R a c A uj c b + u a c A R c b = (17) 
DT a = dT a + uj a d A T d = R a d A e d (18) 

If the torsion vanishes, the latter corresponds to a symmetry property of the Riemann tensor. 
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where J is the generator of rotations. We will see an explicit representation of this algebra 
later in the paper. The left hand side can be interpreted as successive translations along 
e a , — e 6 , — e a , and thus we see that the components of the torsion tensor correspond to the 
non-closure of these successive translations by an extra translation, while the components of 
the Riemann tensor imply that a rotation is also involved. 

In classical general relativity (GR), a basic property of the theory is that the torsion is taken 
to vanish; this is one manifestation of the equivalence principle. In fact, there is a unique 
connection, the Levi-Civita connection & a b, with this property which is determined entirely 
by the co-frame alone (i.e., the metric). Indeed in the familiar Einstein-Hilbert Lagrangian 
formulation of GR, the torsion vanishes as a constraint. In other formulations (the first- 
order or Riemann-Cartan formulations), e a and u a b are regarded as independent degrees 
of freedom and the torsion may then vanish by equations of motion (for suitable choice of 
matter field configurations). In the latter formalism, one can envisage including sources that 
would induce torsion, much as the usual sources induce curvature. It should be emphasized 
though that in our context, we regard e a and u a b as background fields, with no dynamics of 
their own. 

Given the form of the translation algebra (20), the vanishing of torsion in fact corresponds 
to a choice of state. As in the previous section we can consider an elastic medium given by 
a (space-time) lattice A. We will typically be interested in continuum limits, giving rise to a 
continuum quantum field theory, in the presence of a variety of background fields (so that we 
can study various transport properties). At each point in the lattice, we have defined a frame, 
whose magnitudes are tied to the (local) lattice spacing. The commutator of translations 
on the lattice is defined by hopping along a square path; failure to return to the starting 
position corresponds to the path encircling a dislocation of the lattice, and the magnitude 
and direction of the translation determines the Burgers' vector b of the dislocation. There 
exist two primary types of dislocations: (i) an edge dislocation with b perpendicular to the 
tangent vector of the dislocation line (b) a screw dislocation with b parallel to the dislocation 
line (only exists in 3+1-d or higher). An example of the former is shown in Fig. 2b. Now 
consider a continuum limit. If the limit is taken in such a way that a density of dislocations 
b(x) is obtained, we should associate this with non-zero torsion in the continuum theory. 
Lattice dislocations correspond to point sources of torsion. The frame is rotated if the path 
encircles a disclination and continuum limits yielding a density of disclinations corresponds 
to curvature. Disclinations are significant if and only if the field in question carries a non- 
trivial Lorentz representation (that is the generator J ab is non-zero), i.e., it carries spin. The 
effects of dislocations do not carry this requirement. 
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Thus, in condensed matter systems coupled to elastic media, we conclude that the presence 
of curvature and torsion in the continuum limit corresponds to a choice of state. Since both 
curvature and torsion are present, the nature of the background is determined not just by 
the metric, but by both the co-frame and connection. As we will show in detail below, this 
corresponds to the presence of independent Lorentz and diffeomorphism currents (whereas in 
the absence of torsion, these reduce to just the conventional stress-energy tensor). However, 
even in the absence of torsion in the ground state of the system, torsional perturbations 
should be also considered in the context of transport properties. Studying effective actions 4 
of a given field theory in the presence of background co-frame and connections is equivalent 
to studying the correlation functions of the these currents, as the backgrounds correspond 
to sources for the current operators. In some cases, torsion appears in terms in this effective 
action, with coefficients that are physically meaningful. For example, in 2+1 dimensions 
in the presence of time-reversal symmetry breaking, there is a non-dissipative transport 
coefficient (the Hall viscosity) that is the coefficient of a Chern-Simons-like term involving 
torsion, as mentioned in the previous section [6]. 

It is convenient to introduce some additional notation. As indicated above, given a co-frame 
e a , there is a uniquely determined Levi-Civita connection & a b whose torsion vanishes. We 
define the contorsion C a b via 

u a b = & a b + C a b (21) 

so 5 

T a = C a b Ae b (22) 
R\ = k\ + (f)C)\ + C a c AC c b (23) 

Note also that the contorsion is a Lorentz tensor. Generally, we will regard e a and u a b as 
independent. For later use, we will also define 

H = ^H abc e a A e b A e c = e a A T b Vab (24) 

where H abc = —3\C[ a;bc ] = 3T[& c;a ]. H is not in general a closed form, and hence we define the 
Nieh-Yan 4-form 

N = dH = T a A T b rj ab - R ab A e a A e b . (25) 

4 We use the term 'effective action' here interchangeabiy with 'generating functional'. The latter term is 
most appropriate, as indeed, the use of the effective action is that it encodes the correlation functions of 
currents. 

5 We define f) as the LC covariant derivative, (£)C) a b = dC a b + & a c A C c b + C a c A Go c b and k as the LC 
curvature R" \ = dSb a b + Sj a c A & c 6 . 
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4 Generic Properties and Symmetries of Fermions Cou- 
pled to Torsion 



In this section, we will discuss various aspects of fermions on a generic background, mainly 
focussing on the role played by torsion. We are studying Dirac models since they represent 
the minimal continuum models of topological insulators in any dimension. In the following, 
we assume that we are in a d — (D + l)-dimensional space-time with co-frame e a and 
connection u a b with a mostly-plus metric. Results and conventions for spinors and Clifford 
(Dirac) algebra can be found in the appendix. 



4.1 Dirac fermions 



The Dirac action may be written as 6 



1 

~D\ 



e ai...a d e 



ai 



A ... A e D A 



d x det e 



ipj a ip — ipmip 



(26) 
(27) 



We have written the action in this way as it is precisely real (written in other ways, the 
action might be real up to the addition of a boundary term). In odd space-time dimensions, 
m is real, and its sign will play a central role in determining the character of the resulting 
insulating state. In even space-time dimensions m is essentially complex if no additional 
discrete symmetries are imposed (m — > me* 6 * 75 , where 75 is the chirality operator). In 
addition, when torsion is non-zero, there is an additional term 7 that can be added to the 
action, of the form 

S T [e,u) = 4« / deteT a (e b ,eMla,l bc }^- (28) 
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The classical equation of motion for the spinor field involves the Dirac operator 



V = ::{d, + Ajt A + - A u^ bcl bc + B^j + l -aT bc , c - ; 



abc 



(29) 



6 The (Lorentz and gauge) covariant derivative of the Dirac spinor is V^> = dijj + \us a bl ab ^ + Aip, where 
A is an appropriate (non-Abelian) gauge connection. We note also that the invariant form of the action, eq. 
(26), does not involve the frame e a dual to e a . 

7 There are actually two other terms at the same level of power counting. The first, of the form 
i J dete [T a (e i ,e c )ip{'y a ,j bc ]ip — 2Ve a (^"VO] i s Nieh-Yan-Weyl invariant (see below), but a total derivative. 
The second, of the form i J dete T a (e b , e c )ip[ya, J bc }4>, is redundant (it can be absorbed into the definition 
of a £/(!) gauge field). 
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where B a = \T b (e a) e b ) = —\C b a (e b ). The B term arises upon integration by parts in 
deriving the equations of motion. We have included here for completeness a non-Abelian 
gauge field (if the spinor is in a gauge representation t^) and we note that the torsional 
B-term enters in such a way that it looks like it corresponds to an additional gauge field. 
It is not of course independent of the spin connection, but does vanish with the torsion. In 
fact, as explained in [28], the classical theory possesses a corresponding background scaling 
symmetry when m = under which the fields and background transform as 

e a {x) ^ e A(a V(x), uj\(x) h-> u a b {x), (30) 

^ e -(d-l)A(x)/2^ x ^ p ^ e -A( e -(d-l)A/2^ e (d-l)A/2)_ (3^ 

We note that this implies 

T a ^ e A {T a + dA A e a ) (32) 

and hence 

B a = \T\e_ a , e b ) ^ e- A U a + ^e a (A)) (33) 
If we introduce a 1-form B = B a e a , then this is equivalent to 8 

B ^ B + ( ^Y L dA (35) 
which is the transformation of an Abelian (1R + , not U(l)) connection. 

We will refer to this as the Nieh-Yan-Weyl (NYW) symmetry. Note that this is not the 
Weyl symmetry of the metric theory, because in that case, u must transform in order that 
the torsion remain zero. In our case, the Weyl symmetry (at least as far as the Dirac 
operator is concerned) corresponds to a complexification of a U(l) symmetry. In addition, 
the classical Dirac theory also has the usual background diffeomorphism, local Lorentz, and 
gauge symmetries, which we will discuss below. 

Another way to write the Dirac operator is in terms of the Levi-Civita connection, and the 
totally antisymmetric part of the contorsion 

f = (d^ + Ap A + \&^l bc ^j + \c a . tbcl a l bc + B al a + \aT bc]al abc (36) 

= 7X U + Afa + l&^c^ ~ l ~^^ H ^ abC ( 37 ) 

8 Note that the invariance of oj a b implies that the contorsion transforms as 

C\ i \ C\ - dA(e b )e a + V ac i lbd dK{e c )e d (34) 
and the LC connection transforms oppositely. 
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where we have done some 7-matrix algebra (see Appendix) and defined H abc = — 3!C[ o; 6 c ]. 
We will alternately regard H abc as the components of a 3- form H = ^H abc e a A e b A e c or 
as the components of a vector-valued 1-form, H ab = H a b c e c . We note that the parameter 
1 — a determines the coupling of (the antisymmetric part of) torsion to the fermions. Thus 
we can regard it as 'torsional charge', and write 1 — a = q?. For convenience, we will set 
qrp = 1 throughout most of the paper, except in sections 5.3 and 6, where it is illuminating 
to resurrect it. 

Since the Dirac theory is quadratic in fermion fields, the partition function in the quantum 
theory is obtained by performing a path integral over fermions 

Z(A,e a ,u a b ;m) = det(f - m) (38) 

The diffeomorphism, local Lorentz, and gauge symmetries of the Dirac theory remain unaf- 
fected by perturbative (i.e. local) anomalies upon quantization in arbitrary dimension. In 
odd dimensions the NYW symmetry at m = is also non-anomalous. At m 7^ 0, the NYW 
symmetry is explicitly broken. Additionally, the mass term also breaks parity invariance. In 
this paper we will mainly be interested in the quantum effective action for 2+1 dimensional 
Dirac fermions S e ff[e, u, A] = —In det(^> — m). We denote the parity-violating piece of the 
effective action as S oc id[e, lu, A]. In the absence of torsion, symmetry considerations severely 
constrain the form of parity odd terms. For example in d — 3, we have the Chern Simons 
terms 

S odd [e, u),A] = ^ J (cr H A A dA + K H ti (& A d& + ^& A & A &) J (39) 

The coefficient an is called the Hall conductance, while kh, the coefficient of the gravitational- 
Chern-Simons term, has no standard name. Non-zero torsion allows us to construct addi- 
tional terms like 

±J( H e a AT a (40) 

This term was discussed in a slightly different guise in section 2 (see Eq 9); the coefficient 
(h is called the Hall viscosity. Additionally, the effective action also has parity even terms 
of the form 



S even [e, u,A] = ^— J (e abc e a A R bc - ^-H A *H — ^e abc e a A e b A e c ^j 

= ^J d3x det ( e ) ~ : jH abc H abc - 2A^j (41) 



where ^ is the Newton 's constant, A is the cosmological constant, and 7 is a dimensionless 
parameter. We will examine S e ff for the 2+1 Dirac model more closely in section 6. 
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In even dimensions, it is also possible to couple chiral fermions to the frame and connection. 
The action is a straightforward modification of (26, 27) 



1 

D! 



,01 



A ... A e aD A 



1, 



ijj~f ad VP±ip - -V^ ad P ± ip 



d x det e 



(42) 
(43) 



with P± = being the chirality projection operators. The chiral theory also has the 
symmetries of the Dirac theory. However, all the symmetries are spoilt by perturbative 
anomalies upon quantization on generic backgrounds. Later in this paper, we will explore 
such chiral anomalous conservation laws for Lorentz, diffeomorphism, and gauge currents, 
particularly in 1 + 1 dimensions, and their connection with S Q dd for the 2+1 Dirac model. 
We will see that while torsional terms like (40) leave consistent anomalies unaffected, they 
do modify the covariant anomalies. 



4.2 Classical Ward identities 

In this section, we state the classical conservation laws for fermions coupled to the coframe, 
connection and a U(l) gauge field. 9 Although we will discuss these in the context of 
Dirac fermions (for arbitrary d), the results generalize in a straightforward manner to chiral 
fermions in even dimensions. Let us begin by defining the following currents 

(J/*) = g^e^ (44) 
W = ^7 a V^-V^ 7 >) (45) 

(JT b = \f£H* h 1> (46) 

which couple respectively to the gauge field, coframe, and connection in the classical action. 
In the absence of torsion the last two currents are not independent. The components of the 
current J a give the usual notion of the stress-energy tensor via 

T,u = J a ,e b u Va b (47) 

Also note that the spin current J^ b vanishes in d = 2. It will be convenient to introduce the 
corresponding 1-forms J = J^dx 11 , J a = J^dx^ and J ab = J^dx^. Invariance under U(l) 
gauge transformations implies that J is conserved, i.e. d * J = 0, which in components is 
the usual d fJi (det(e) J M ) = 0. 

9 In the rest of the paper, we will restrict the gauge group to U(l) in favor of somewhat simpler notation. 
We will use the symbol q for the Z7(l) charge. 
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4.2.1 Diffeomorphisms 



The invariance of the classical action under local background diffeomorphisms follows im- 
mediately from writing it as the integral of a top form, as in (26). We will take the action 
of local diffeomorphisms on fermions and background fields as 

<^ = i e W, Se a = D£ a + i{T a , 5uj ab = i^R ab , 6A = i^F (48) 

where £ is a vector field with compact support and i% is the interior product of £ with a 
differential form. These transformations differ from ordinary diffeomorphisms by local gauge 
transformations, so we will refer to these as covariant diffeomorphisms. Using equations of 
motion for the fermions, the variation in the action under (48) is given by 

5Diff.S = J [i^FA*J+(D^ a + t^T a )A*J a + t !; R ab A*J ab ] (49) 

and so invariance of the action implies the classical Ward identity 

D*J a - ieaT b A *J b - i^Rbc A * J bc - i*F A * J = (50) 

4.2.2 Local Lorentz transformations 

The spinors and background fields transform under an infinitesimal Lorentz transformation 

as 

# = Se a = -e\e\ 5u\ = -{D6)\ (51) 



Under (51), the action changes by 

hor.S = - J [D6 ab A * J ab + 6 ab e a A *J b ] 

The Ward identity is 



(52) 

D*J ab - e [a A * J b] = (53) 



4.2.3 Nieh-Yan-Weyl transformations 



The action on fermions and background fields is given by 



5i) = -- — Se a = Ae a , Soj ab = (54) 



Under Sip = —^^Aij;, the action transforms as 



5S NY w = -(d- !) y A [Vabe a A*J b -m vol ipip] (55) 
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The second term, where vol is the volume form, is present because the mass term explicitly 
violates the NYW symmetry. For m = 0, we have the Ward identity 

Vab e a A *J b = (56) 

In components, this is T M M = g^e^Jfyab, the trace of the stress-energy tensor. Thus in this 
sense, the NYW symmetry gives rise to the same conservation law as does Weyl invariance 
of the second-order formalism. 

4.3 Lichnerowicz-Weitzenbock Formula 

The heat kernel for the operator Tp 2 will play a central role in some of our computations, so 
before moving on, let us briefly describe this operator. We begin by noting 

P 2 = -f(y a + B a )^(y b + B b ) (57) 
= 1 a 1 b (D a + B a )(D b + B b ) (58) 
= ri b V a V b (59) 

where D a is fully (Lorentz) covariant and T> a = De a + B a . In manipulating this expression 
we need various facts about the Clifford algebra (see Appendix) and we also encounter the 
commutators 

[D a ,D b ] = -T: b D c + ^R cd , ab7 cd + tqF ab (60) 

[D[ a , B b] ] = - l -T: b B c + l -G ab (61) 

where G^ u = d^B u — d u B fl . Consequently, the Lichnerowicz-Weitzenbock formula takes the 
general form 

f 2 = rj ab V a V b -\r+ ^F abl ab + \R cd ; ab l ahcd + \l ab G ab - \i ah Tl b V c - l -R b a . <dbl i&2) 

In the absence of torsion, the curvature tensors satisfy R ab ;cd = Rcd-, ab and R b a . bd = R b d . ba . 
Therefore the last four terms in (62) would vanish in the torsionless case. 

5 The Hall Viscosity 

Here, we will give an extensive discussion of the Hall viscosity for the Dirac model in 2+1. 
We begin with 3 separate derivations (which are of course equivalent): the Berry curvature 
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due to area-preserving diffeomorphisms, linear response via the stress-stress correlator, and 
the effective action in background fields. In each of these three subsections, we will obtain 
the same result, written in terms of an integral It(ui), which contains a power divergence 
at high-energies. In the follow-up section, we describe how to interpret this divergence in 
terms of a Pauli-Villars regularization scheme. 



5.1 The Berry Curvature 

We begin with the direct calculation of the Berry curvature for a Dirac field on a spatial 
torus with complex structure r. This calculation is the analogue of that done by Avron et 
al. for IQHE[4, 23, 8, 29]. 

We consider a square torus, made in M 2 by identifications (x, y) ~ (x + aL, y + bL) with 
a, b G Z. We take this to have fixed volume L 2 , and consider area preserving diffeomorphisms, 
which we take to correspond to modifications of the metric of the form 



« = 4G kV>' = ^ (63) 



The basis vectors are then 



T2d x , e 2 = -— (-nd x + d y ) (64) 

T 2 



e 1 = ——{dx — Tidy), e 2 = y/r^dy (65) 



and the veirbein is 

e x = - 
/r 2 

Since the components of the metric are constant, we will take the connection to vanish. Take 
a basis for C£(2, 1), 7 = ia 3 , 7 1 = o\ and -y 2 = a 2 . The Dirac operator is 

tid -t= (fd x — d v )\ 

and the Hamiltonian is then 

' m p 



H = _ H (67) 
yp —ml 

where p = ^= (fp x — p y ). The eigenstate energies thus satisfy 

P 2 2 || 112, 2 \ T \ 2 Pl -2Tip xPy +p 2 

ti = pp + m = \\p\\ +m = + m . (68) 

T"2 
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We introduce the phase rf = = y |. States in the negative energy band can then be 
written in the form 



lMfc,*,;r)=| y^L_| 



-77 



2\E\ 



The states in the upper energy band are 



\E\+m 



Mp*,Pv;t)=\ , v ^| (70) 



^ V 2\E\ 

with E given by (68). On the torus, the components of momentum are of the form 

p x = 2n(q + h/2), p y = 27r(r + k/2), q,reZ (71) 

where (h, k) label the spin structures. 10 We will consider the insulating ground state in which 
all of the negative energy states are occupied and the positive energy states are empty: 

i n > = II4wi ± »s ^ 

k 

where V^{r) is the projection operator onto the negative energy state at each momentum k 
for modular parameter r. From the above considerations, we find 



The Berry connection is 



' \E\-m _ l \E\+m 

v % {f) = ( '' V v y^m ) (73) 



a = -i(n\d\n) (74) 

= -^J2 trV l^ dv k( r ) ( 75 ) 

k 

= ^i(q,r;r)dtp-(q,r;T) (76) 



In writing this expression, we have explicitly assumed that the oscillators have no r-dependence 
and are normalized via {dt r , d q > jr >} = ^5 q -q>6 r - r i. A straightforward calculation then gives 

A = ^J2f(M 2 M) 2 d(v 2 ) (77) 

q,r& 



3 That is we require 

i/>(x + a,y + b)= e lw( - ha+kb) ij(x, y), a,beZ. 
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f(\\p\\ 2 ) = m (7 8 ) 



where 

f(\M\ 2 \ = - 

I 2 

We find the Berry curvature from taking the exterior derivative 

i E f(\\p\\ 2 )dpAdp (79) 



% 

m,n(zZ 

dr Adf 



= i- 



2t l 

q,r 



£*£/'(lbll 2 ) ( 8 °) 



This discrete sum can be cast as an integral in the large volume limit 

£ = E^'(INI 2 )^ 2 / (^AlbllVx (si) 



We introduce polar coordinates by transforming 



V^Px = ||p||cos0, —^(-Tip x +p v ) = \\p\\smO (82) 

V T 2 



so the integral becomes 



S = (2^/ dpp3f '( p2) ^ d9cos29 ( 83 ) 



L 



2 



dy yf(y) (84) 



87rr 2 jo 

L 2 r , d m 



dy y- (85) 

= 8^ /T(m) (86) 

This defines the integral Irim) which diverges at short distances (large ||p||). Thus, the 
calculation should be performed with an explicit regulator, a subject we will address care- 
fully below. In the calculation performed by Avron et al. for the IQHE produced from 
Landau levels the result is finite and proportional to 4- where £b is the magnetic length; 
a short-distance cut-off generated by the uniform magnetic field. For the Dirac insulator, 
and topological insulators in general there is no such natural length scale which is why the 
calculation must be carefully regularized. 



5.2 The Stress-Stress Correlator 

Now, let us consider the 2+1 Dirac theory and look directly at {T^ u T\ p ). By general argu- 
ments, the Berry curvature computed above can be related to (the parity odd part of) this 
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correlator, and in particular to the Hall viscosity. We note that a careful analysis of the 
Kubo- formula description of Hall viscosity transport, and its connection to conductivity for 
Galilean invariant systems, is discussed in Ref. [11]. In any case, we show here the explicit 
calculation for clarity. Recall first the Lagrangian 



m 



i=x,y a=l,2 



i7) 



The stress tensor is then given by 



T 0i (t,x) = - : $(t,x) [ida°] i/i(t,x) : 



we compute the correlator 

(T^frf^OjO)) = -Atr 



iSp(t; x) 



?9) 



iS F (-t;-x) (90) 



where iSf = (i</D — m)^ 1 . The 8ij term in 7^ will not contribute to what we are interested in 
(although it would contribute to other visco-elastic response coefficients). Since we are just 
interested in this response, we now set = 5?. Now we can Fourier transform 



(TijT ki )(p) 



-Atr 



d 3 Q [|?(i7j)] d + m ) [\ r {kli)] (f + m) 



(27r) 3 (q 2 — m 2 (q 2 ))(r 2 — m 2 ) 
where f=p + q. The terms that have an Levi-Civita tensor are the following 

d 3 q 



(91) 



rn 



(2ir) 3 (q 2 — m 2 )(r 2 — m 2 ) 
d 3 q m 



q{ir(k tr 



(27r) 3 (q 2 — m 2 )(r 2 — m 2 



(92) 
(93) 



This can be rewritten 
d 3 q 



^7° + Pm7 m , q 2 = - <f ) 



2ie 



Mi 



mqi)r k )U q 



(27r) 3 (q 2 — m 2 ){r 2 — m 2 ) 



2ie 



(t(3 



d 3 q 



mqi)r k )U r 



(27r) 3 (q 2 — m 2 ){r 2 — m 2 ) 



(94) 



where we used tr^jijj = —2ieij in all terms. Now we expand this to zeroth order in the 
external momentum p since it is already first order in the external frequency, and get 



2imuj n e 



d 3 q 



Qi)Qk) 



(2ir) 3 (q 2 — m 2 )(r 2 — m 2 ) 



(95) 
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Here we used (cj q — u r ) = — cu p . We can do the u q integral, which is trigonometric: 

1 f d 2 q me {eu q i)qk) 

2 P J (27r) 2 (<f + m 2 ) 3 / 2 1 J 
1 f d 2 q _g d m 



4 p ^ w 7 (2tt) 2 " <9,f (^ + m 2 )V2 
= -T^u P £{£{j$i)k)lT(m) (97) 

l07T 

This is the same integral that appears in the Berry phase calculation. One can also similarly 
establish 

= ^PnS i(j e k)n J o dyy^ (y+ ™? (y))1/2 (98) 
and thus for small we have 

(T pu T Xp )(p) = ■^5 {u{p e^ X )aP a lT(m) (99) 



5.3 The Background Field Formalism 

This result in Eq. 99 can also be rewritten in terms of the stress current (as above, a 
vector-valued 1-form current associated with diffeomorphisms) 

(J a A)iP) = 7^ W^Mrn) + • • • (100) 

P l07T 

At leading order, this may be interpreted as coming from a quadratic parity-violating term 
in the effective background-field action of the form 

S odd [e,co,A}= ^... + J_/ T ( m ,A) J e a Ade b Vab + ..}j (101) 

In addition to this term, there will be a higher order term that renders 11 it locally Lorentz 
invariant 

S odd [e, u, A] =(...+ ^/ r (m) J e a A T b Vab + (102) 

This result means that the Dirac fermion contributes ^-Ix(m) to the Hall viscosity (h- 
Again, we stress that the result is divergent, but we will consider this divergence carefully 
later in the context of a Pauli-Villars regulator. 



11 In a generic state of a condensed matter system, we of course have no reason to suppose that Lorentz in- 
variance would be obtained. However, the specific model under discussion does have local Lorentz covariance, 
and thus the effective action will share that feature. 
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It is instructive to produce (102) directly via a background field expansion. Here what we 
will do is consider the Dirac fermion in a general background co-frame e a and connection 
u a b and integrate out the fermion. This is in principle straightforward, since the fermion 
integral is Gaussian and the fermion is gapped, although the result is non-standard, as we 
will consider a torsional connection. Later we will extend the calculation to include the 
effects of curvature, as well as a calculation of the parity-even terms. 

We have written the general Dirac action in a previous section. Here we are interested in 
computing the determinant of the Dirac operator (now for a general value of q?) 

f = l ae a (®n "I" A^t^ + ~^;fec7 bc ^ + Bf,"y b + ^<yT a . bc 'y abc (103) 

= Tp + \q T T a , bcl abc (104) 



Irpb 

2 1 al 

Formally, we write this as 



where B a = \T b b . In particular, the contribution to the effective action is — lndet(^) — m). 



o 77) _ m 

-lndet($> -m) = - \ndet(f - m) - lndet \ (105) 

lp — vn 

o 

= - In det(f-m) -lndet \ £ (106) 

Tf> — mTp + m 

= - Indetpp - m) - lndet ^ " \f ~ + (107) 

Tp — m 2 

o 

where by Tp, we mean the Dirac operator with (torsionless) Levi-Civita connection u , while 
Tp is a general torsional connection. The second expression vanishes in the absence of torsion 

o 

while lndet(^ — m) is the usual gravitational effective action (which will contain a volume 
divergence, the gravitational and gauge Chern Simons terms, as well as higher order terms). 
See Ref. [30] for related calculations. 

Let us then focus on the second expression. One advantage to this quantity is that the 
volume divergence (being independent of the connection) has to cancel exactly, because 
the numerator and denominator differ only in the connection (and not the frame). The 
expression is well-defined mathematically, as it is comparing the contributions of two different 
connections (analogue of relative Chern-Simons). 

We have (as a special case of (62)) 



TP 2 = V ab VaV b -\°R+^F ab r b (108) 
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and one can establish 12 



f-f = \c bc {e a )l a l hc + B bl b + l -aT a , bcl abc (109) 
= f C bc (e a h abc = q -fe abc C a , bc ee \q T c (110) 



Formally, we may write 



-lndetf — — = -lndeth+fc^-m)- 1 ) (111) 

~ -^-tr c(f -m)' 1 + ... (112) 

So the entire effect of torsion can be thought of in terms of a cipip vertex, and the above 
expression just corresponds to the fermion loop diagram with a single insertion of c (the 
singlet part of the contorsion). It will be sufficient to work to first order in c to see the 
torsion Chern- Simons term, as 

e a A C ab A e b = e a A T a = d 3 x e^e^C^e^ = -cd 3 x (113) 

where in the last equality we expanded around flat space-time. Indeed, to extract just the 
torsional CS term, we can expand around flat space-time to obtain the leading result 

- \m Tr c(D - vr?Y 1 = \m J ' d*x c(x) j ^ f - i - (114) 

= Urn [ d 3 xc(x) ds [ ^L e -^+mi)s (n5) 
2 J J e J (27r) d 

= ^m/(0,0; \m\) J d 3 xc(x) (116) 

Here we are using the notation 

I(p, g] \M\) = j±j £^e-^ 2+M ^\k\« (117) 

Now we can do a Wick rotation in the complex k plane (k = —ik^) and then rotate the 
contour in the complex s-plane (s = —it) to get 

I(p,q; \M\) = -^l§j ^e"^ 2 )^ 2 (118) 

12 Here we have a choice for the representation of the Clifford algebra, namely 7 abc = ±e abc , which essen- 
tially amounts to a choice of orientation. We will choose the positive sign; this will be reflected in the sign 
of (h , and later on other parity odd coefficients. 



28 



We note that the A;- integral converges in the UV unless t = 0. So we cut off the t-integration 
at t = e ~ f/A 2 . 

Iip^M) = -^{^)\M\ 2 ^ +1 t{- P -\-\aM\^ (119) 

In terms of the previous notation, one finds ir(wi) = 87rm/(0, 0; For later use, we note 
that for small e, we have 

zl(-l,0;m) = --J— + ... (120) 
87r|m| 

/(0 ' 0;m) = -T,vm + i^ + - < 121 > 

2 1 \m\ 2 1 Iml 3 

' i/(1 - 0im > = 3(i^-^vs; + ^r + - (122) 

There is also a Chern-Simons term for the U(l) gauge field, including which we get 

S odd [e, u, A] = ~q T mI(0, 0; |m|) J e a A T a + q 2 m/(-l, 0; \m\) J A A cL4 (123) 

Comparing these results with (39) and (40), we get the unregulated coefficients for a single 
Dirac fermion of mass m and charge q 

( H ( m ) = -q T mI(0,Q;\m\) (124) 
a H (m) = -2iq 2 ml (-1,0; \m\) (125) 



5.4 Pauli-Villars Regulator 

Now, we will finally consider a proper regularization for the Hall conductivity and Hall vis- 
cosity. To do this we introduce a set of A" regulator fields with masses Mi and normalizations 
Cj. The regulated coefficients are then 

N 



( H ({M t }) = -9*^6^1(0,0; |Mi|) (126) 

i=Q 
V 

a^({Af ?; }) = -2^ 2 ^aAf,/(-l,0; 1^1) (127) 
where the original fermion is labeled with Mq = m and Cq = 1. 
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What we should do here is choose the remaining Cj and Mj such that the resulting expressions 
are finite as e — > 0. Making use of the expansions (120-121), we find as e — > 



Cff({M 4 }) 



A? 



AT 



+ <7 i |M| 2 + 



7T6 



47T 



(128) 
(129) 



i=0 



where <7j = sign(Mi), and in particular cr = sign(m). To render these expressions finite, we 
must require 



N 



J]aM 



0. 



(130) 



i=0 



We note that in fact this condition must always be satisfied regardless of whether or not the 
system is a trivial or a topologic insulator, that is, it is satisfied independent of the sign of 
m. To see that this can be achieved, we take as a guide the (2+l)-d lattice Dirac model [12 
with Hamiltonian 

H = ^2 = ct a [(m + fi bw (cos k x + cos k y - 2))er 3 + v F sin k x a 2 + v F sin fe^eri]^ c^. b (131) 

We compare these continuum calculations to the lattice Dirac model which can be interpreted 
as having one low-energy fermion and three regulator fermions. On the lattice, momenta 
are measured in units of the lattice spacing a, and the Brillouin zone is a torus given by the 
square [— f , ?] x [— - , ^1 with the boundaries identified. The spectrum can be interpreted 
as four independent copies of the (2+l)-d continuum Dirac Hamiltonian located at four 
separate lattice momenta k = (0, 0), (0, -),(-, 0), (-, -). The Dirac fermions away from the 
origin can be thought of as analogous to the Pauli-Villars regulator fields, and the regulated 
continuum model is then an approximation to the lattice model where only the effects of 
these Dirac points are taken into account. Examination of the Hamiltonian near these points 
gives 



Mi 


Q 




m 


+ 


m 


2/^(m) 




m 






m 




+ 



and we note that indeed ^2^ =0 CiMi = 0, independent of the value of 



m. 
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Returning to the general Pauli-Villars problem then, we will take N = 3 with 13 Cj = 
{+, — , — , +} and then take masses such that M1 + M2 — M3 = m, say by taking Mi = m+Ai, 
M2 = m + A2, M3 = m + Ai + A2, where we are considering m << Ai^. Now, if we use 
<jh as a guide, we believe that the sign of m determines two phases, one of which is a trivial 
insulator and one of which is a non-trivial topological insulator. A choice of the signs of 
A 12 will correspond to choosing which value of cr = sign(m) gives the trivial phase. For 
definiteness then, let us choose Ai 2 > 0. We thus obtain 



Ch ■ 
a H ■■ 

Thus we see that 



Q.T r 
47T 



[o"o| 



m 



(m + Ai) 2 - (m + A 2 ) 2 



m 



+ Ax + A 2 ) 2 ] 



m 



+ 2A!A 2 ] 



Air 



K - 1 - 1 + 1] 



(132) 
(133) 



m > 0: 



m < 0: 



Ch 
oh 

Ch 
oh 







2Z Al A 2 + 2Z m 2 

Z7T Z7T 

2?r 



(134) 



The value of Ch diverges as we remove the regulator masses to infinity, but in a way that 
is independent of the phase. Thus, including a counterterm to subtract this divergence 
(or by including additional regulator fields appropriately), we see that (the renormalized) 
coefficients Ch and an vanish in the trivial phase, while in the non-trivial phase, they are 



q T 



m 
2^ 



(135) 



Said another way, it is the difference of C,h (as well as an) in the two phases that is universal. 
In fact, this difference can be detected by studying anomalous conservation laws for 1+1 
dimensional chiral fermions which live on the edge separating the trivial and non-trivial 
phases, as we will see shortly. We give a more physical discussion of the Hall viscosity in 
terms of the bulk-boundary correspondence in Section 9. 



13 Since each of the Dirac points has a Hamiltonian linear in momentum, Cj refers to the "parity" of the 
Dirac copy; thus the modes at (0, it), (jr, 0) have opposite parity to the other two. Parity can be understood 
as the relative sign between the two terms in the Dirac Hamiltonian which are linear in momentum. In the 
following section, we will see that the condition X^ilo — w iU follow from finiteness in the parity even 
sector. This is the usual "fermion doubling" phenomenon. 
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6 Effective action as Chiral gravity 



In the previous section, we have shown how to extract and regulate the Hall viscosity term, 
working in a flat background. In this section, we compute the effective action for massive 
Dirac fermions coupled to a generic background frame and connection. We will work in the 
limit of large radii of curvature and torsion in comparison with the fermion correlation length 
£ ~ — . We will utilize the asymptotic expansion for the trace of the heat kernel of Tp 2 (see 
Appendix A) to extract terms which survive as (ma) — > oo (a being the radius of curvature). 
The resulting renormalized action in the non-trivial phase differs from the trivial phase by 
an SX(2,R) Chern-Simons term. 

We will work on a manifold M with Euclidean signature. At this point, we might expect 
that the action can be written in the form 

1 f (ia H A A dA + ik h CS[u] + i( H e a A T a + —e abc e a A R bc - —vol M + ■■) (136) 

2 J \ - K N K N J 

where CS[u;] = uj a b A duj ha + \oj a b A u b c A u c a and voIm = ^abc^ a A e b A e c . More precisely, 
we will eventually find that the leading terms in the (renormalized) action will organize into 
a specific form involving e a , u ab and H. 



6.1 Parity odd terms 

In the absence of background torsion, the parity odd effective action for the Dirac theory is 
given by [30, 31] 14 

S odd [e,uj,A) = ml (-1,0; \m\)q 2 J A A dA + ^ml (-1,0; \m\) J tr(# A d& + ^& A& A&) + ... 

(137) 

where the integral I(p, q; \m\) has been defined in (117). In order to extract terms to leading 
order in torsion, the simplest thing we can do is repeat the relative Chern-Simons calculation 
of the previous section in the presence of background curvature. We thus replace (114) with 

i ° 2 _ 

S odd [e, u, A] = S odd [e, u, A] + -q T m Tr c(-Tp + m 2 ) 1 + ... 

= S odd [e, Cj, A] + -q T m / dt Tr ce"'^ W) + ... (138) 
4 Jo 

o 

14 We will assume for convenience that Tp and lf> have no zero modes. 
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where as before, c = e abc C a -b c . In the limit t — > 0, there exists an asymptotic expansion for 
Tr e*^ (in arbitrary dimension d) 

^ e *^(4^I> ft (139) 

where are integrals over M of polynomials in curvature, torsion and their covariant deriva- 
tives. The important point is that it suffices to use this asymptotic expansion in order to 
extract terms which survive as (ma) — > oo. For the case at hand, we have (see Appendix A 
for details) 

Tr / * J (l-±R + 0(t 2 )) vol M (140) 

where voIm is the volume form on M. The 0(t 2 ) terms are unimportant for the present 
discussion, as they are higher order in powers of curvature and lead to negative powers of 
(ma). From (138) and (140), we get 

S odd [e,u,A] = - l -q T ml (0, 0; \m\) J e a A T a + ml (-1,0; |m|)g 2 / A A dA (141) 



+ — m/(-l, 0; \m\) J [ti (& A d& + -& A & A &) - q T R e a A T a + ■ ■ 

where the ellipsis indicates terms higher order in torsion. Since we are considering this result 
as a functional of u, we can re- write it in those terms 

S odd [e,oj,A\ = - % - q T mI (0, 0; \m\) J e a A T a + ml (-1,0; \m\)q 2 J A A dA 

+ ^m/(-l,0; \m\) J ^tr (u A du + A u A u)j H (142) 

up to terms quadratic in curvature and torsion. 

It will also be convenient to introduce a new 5*0(3) connection 0J^l b = ti^-ob ~ f-H^;o6- This 
is natural because, as we have seen, the Dirac operator involves only u and H. If we choose 
the specific value /3 = —qr, one finds that the second line of (141) combines into a single 
Chern-Simons term for oj^.^ ■ In fact, this is indicated by the structure of the chiral anomaly 
in d = 4 (see Eq (A. 7) in Appendix A). 15 Given the level of the calculation that we have 
presented, it is possible to verify this to linear order in torsion. Equation (141) may be 
rewritten in these terms (to linear order in torsion) as 

S odd [e,u,A] = - l -q T mI(0,0;\m\) j ' e a AT a + q 2 ml (-1,0; \m\) J ' A A dA (143) 

+ J-mI(-l, 0; |m|) I tr (u/~ 9t) A dtu ( - qT) + -w ( " 9t) A u^ qT) A co ( - qT) ) 
24 / 3 



15 The connection between the d = 2n+ 1 parity odd effective action and d = 2n chiral anomaly is provided 
by the Atiyah-Patodi-Singer index theorem; see [31] for details. 
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6.2 Parity even terms 



The parity even terms in the effective action are given by 



f°° dt 

S even \e, U), A] = lim / — Tr 



-tm 2 +fp 2 



(144) 



Once again, it suffices to use the asymptotic expansion for Tr e*^ in order to compute terms 
which survive in the large (ma) limit. The asymptotic expansion in this case is given by (see 
Appendix A) 

Tr * / (l - + (t 2 )) vol M (145) 

where R^~ qT ^ = R — ^-H a b c H abc is the scalar curvature constructed out of uj^ . Using 
(144) and (145), we get 



S, 



r(-qt) 

lim / (iJ e (l,0;H)+/ e (0,0;|m|)— — + -•• 



12 



vol 



M 



(146) 



= hm y"^J £ (l,0;|m|)t; O / M -^/ £ (0,0;|m|)e abc e a Ai?(-^ c + --^ (147) 
where the ellipsis indicates terms of order (ma) -1 . 

In order to regulate the divergences in the effective action, we introduce N Pauli-Villars 
regulator fermions as we have done previously. Using the expansions (120-122), we find 



N 



i=0 



Tie 



+ cr l \M i \ 2 + 



N 



a H ({M l }) = -J2C l a l + --- 

i=0 

1 - 

({M,}) = — £cta + ... 

907T L — 4 



i=0 



1 1 N 



2k n 



i=0 



A 



iV 



-(W) = £a 

1=0 



7T6 



|Mt| + • • 

|M f | 

+ 



+ 



_3(4vre) 3 / 2 4vr0F^ 6tt 
We require that the terms that diverge as e — > have zero coefficients. This implies 

N N N 



(148) 

(149) 
(150) 
(151) 
(152) 



£q = o, £(7^ = 0, E^I Mi 







(153) 



i=0 



i=0 



i=0 
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Thus, we have one new condition from the parity even sector and we now see that the first 
condition (we used this above) is also required by the parity even sector. If we assume for 
simplicity that all of the regulator masses are positive, 16 then we arrive at 

Ch = VT^^- (154) 

°H = -73 7T- (155) 

Z7T 2 

KH = (156) 

J_ = H 1^0 {r) 

k n 12tt 2 y J 

A N 

— = 4ew ( l58 ) 

k n 6vr^ 

where again, <jq = sign(m). If we examine the conditions (153), we can furthermore establish 
that 

— = Al - - m 3 — — ^ (159) 

KAT 37T Z 

for a quantity Ao that generally scales with the regulator masses, but is independent of 
a . We thus see that apart from the Aq term, all of these coefficients vanish in the trivial 
phase (cr = 1) and the effective action there is just S^j = Aq f vol M , a pure cosmological 
term (of course, there are also higher order terms in curvature and torsion, which decay 
exponentially or as negative powers of (ma), that we have not included here; those terms 
then determine the transport properties of the trivial phase). The non-trivial phase has an 
action consisting of the same Aq volume term, plus an action that is known as chiral gravity 
(as well as the usual U(l) gauge Chern-Simons term). In other words, the difference of the 
gravitational actions between the two phases can be written in terms of the Chern-Simons 
form of a single SL(2, K) connection. 17 Indeed, writing u a = ^e a b c a/ c , we define a connection 
j[ a = w (^)' a — ije a . One then finds 

iCS[A a ] = 2i (^A a A dA a - ^e abc A a A A b A A c ^j (160) 

= —%CS[iO^\ + ^ v °l + j£ a bce a A R^' bc — i^e a A (161) 
= -iCS[J^] + ^vol + % c e a AR^ bc -Me a AT a (162) 



16 This assumption leads to the m > phase being trivial. Another choice would make the m < phase 
trivial. 

17 Here we are using the language of real time. The group theory involved here is that the isometry group 
of AdS 3 is ~ SO(2,2) ~ 50(2,1) x 50(2,1) ~ SL(2,R) x SL(2,R). See Ref. [3 ] for details. 



35 



and we thus see that if we identify = —qx as above, £ = (2\m\) 1 and A = — 1/£ 2 , the 
action in the non-trivial topological phase is 

S; ff = St ff + ^CS[A*] - ^CS[A] (163) 

The Chern-Simons level k evaluates to 1/24. 18 Incidentally, chiral gravity has been studied 
in the context of holography [3 3], in which the gravitational fields are dynamical. Indeed 
if we introduce the notation /i = 2 kI,k h ( nere = — 1); the Brown- Henneaux formula in 
asymptotically- AdS^ geometries gives the central charges of the dual 1 + 1- dimensional theory 

as 

(, l \ Vivify 1\ MnA . 
cl = (l- — ),c R = 1 + — . (164) 



Thus in the holographic case, we have Cf, = 1, cr = 0. In the present case, this 1 + 1- 
dimensional matter is supported on the interface between the topological insulator phase 
and the trivial phase. In the following sections, we will investigate this further and show 
that the parity odd transport coefficients of (163) are encoded in anomalies of the 1 + 1- 
dimensional theory on the interface. This will include not only the chiral anomaly, relevant 
to the charge sector, but also diffeomorphism and Lorentz anomalies for the gravitational 
sector. We will see that indeed (163) implies that cl — cr — 1. 



7 The Callan-Harvey Anomaly inflow 

Consider the non-trivial phase labelled by non- vanishing parity odd coefficients (cr#, (h,kh) 
on a 2+1 dimensional manifold M, separated from the trivial phase by the 1+1 dimensional 
boundary S = dM. This can be thought of in terms of a 2+1 dimensional Dirac fermion 
with mass m < on M, and m > outside, with some interpolation region, the interface S, 
which we refer to as the domain wall. In general, there could be multiple fermions with mass 
domain walls along S, and their number decides (cth,(b, k h)- The domain wall hosts 1+1 
dimensional chiral fermions, whose anomalies will encode the shifts in (an, Ch-,^h) between 
opposite sides of the domain wall[ ]. In the absence of curvature (we will return to the 
general case later), the parity odd effective action can be taken to be 19,20 

S oMM i k [e, oj, A] = ^ / e A A T A + ^ / A A dA (165) 

18 This result satisfies the quantization condition k G jgZ given in [ ] for manifolds which admit a spin- 
structure. Here we get twice that result, because we have a full Dirac fermion in 3d. 

19 The gravitational Chern Simons terms (proportional to kh) lead to currents which are proportional to 
the Levi-Civita scalar curvature. Hence we ignore these terms temporarily. 

20 In this section, we will use upper case letters for Lorentz indices in the bulk and lower case letters for 
Lorentz indices on the boundary /domain- wall E. 
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Let us first focus on the gauge Chern-Simons term and review its relationship with anomalies 
in the boundary. In the presence of a boundary, the U(l) Chern-Simons term is diffeomor- 
phism and Lorentz invariant, but not gauge invariant. Under a gauge transformation we 
have 

S a S dd,buik = / da A F = / aF (166) 

Gauge invariance implies that this should be accounted for by the U(l) anomaly of chiral 
fermions localized on S. For left-handed and right-handed chiral fermions on the 
edge, the anomaly is given by 

6 a Sx = / aF (167) 

This cancels the variation of the bulk action provided q 2 (riL — tlr) = —2tt<7h, which is 
indeed the case as can be checked by constructing the localized zero modes of the bulk Dirac 
operator (see [1] for details). The anomaly in (167) is called a consistent anomaly, because 
it is obtained by the variation of the chiral effective action in 1+1 dimensions. We refer to 
the corresponding non-conserved current as J cons , with the anomalous Ward identity 

a * Jcons = — -: q F = — F (168) 

Returning to the bulk, the variation of the effective action with respect to the gauge field 
determines the current 

SS oddMlk = a H f 5A A F + f 5A A A (169) 
We can read off the bulk U(l) current from here 

* Jbuik = oh F (170) 

which is conserved by virtue of the Bianchi identity, i.e. d * Jbuik = 0. However, the flux of 
the bulk current into S is non-trivial and is given by 

AQ S = J *J bulk = a H J F= J F (171) 

We can interpret this as the charge injected into the edge from the bulk, but notice that it 
is twice as much as the consistent anomaly in (167). To explain this apparent discrepancy, 
notice from (169) that there is an additional boundary current induced from the bulk 

*3 = °fA (172) 

This prompts us to define the net boundary current J cov = J CO ns + j, which we will call the 
covariant current. The conservation equation for J cov is now 

d* J cov = a H F= T ^^F (173) 

Z7T 
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which agrees with (171). The anomaly in the form (173) is called the covariant anomaly. 21 
We see therefore, that the covariant current in the boundary carries the charge which is 
injected into it from the bulk. 

The Hall viscosity term in (165) on the other hand is invariantly defined under diffeomor- 
phisms, Lorentz, and U(l) gauge transformations, and thus does not lead to consistent 
anomalies in the edge theory. The consistent frame current J^, ns in the edge is therefore 
symmetric, and suffers only from the anomaly due to the U(l) Chern-Simons term 

D * J a cons - ie«T b A * J b cons - i e »F A * J cons = - q 2 ^-^( leaA )F (174) 

e [a A * Jl s = (175) 

(Recall that since the domain wall is 1 + 1-dimensional, the Lorentz current J ab vanishes.) 
Equation (174) is not gauge covariant. However, it is clear what we must do - we shift to 
the covariant currents. 

Consider then, the variation of the bulk action under a change in the frame and connection 
SS oddMlk = ( H [ 5e A A T A - % f 5u AB e A A e B + % f 5e A A e A (176) 

JM * JM * JT, 

In the boundary term, we should interpret the result in terms of fields defined on the bound- 
ary. Generally (as was implied in the discussion of the gauge case above), p-forms will pull 
back to the boundary. In the case of vector- valued forms e A and u AB , we also must decom- 
pose the pullbacks in representations of the boundary Lorentz group. Generally, we are free 
to choose independently a co-frame E a and spin connection Q ab in the boundary. These can 
be identified with the pull-backs of e a and u ab up to a Lorentz transformation. The normal 
components e n and u na represent extrinsic effects. Conventionally, the pullback of e n to E 
vanishes, which can be achieved by a local bulk Lorentz transformation of the frame. 

We read off the bulk frame current and spin current from (176) 

*JL, = ChT a , *J^ k = - C fe A Ae B (177) 

21 The reason for the terminology consistent and covariant comes from the more general case of non-Abelian 
gauge anomalies. In that case, the consistent anomaly satisfies the Wess-Zumino consistency condition, but 
fails to be gauge covariant (it involves dA rather than F). The covariant anomaly on the other hand, does 
not satisfy the Wess-Zumino consistency condition, but is fully gauge covariant. The consistent and covariant 
versions of the anomaly differ by current redefinitions which do not come from local counterterms, but are 
equivalent so far as anomaly cancellation is concerned. The difference between the covariant and consistent 
currents is usually referred to as the Bardeen-Zumino polynomial. 
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while the frame current induced in the edge theory is 

*r = Y ea ( 178 ) 

It is easy to check that the bulk currents satisfy the proper (non-anomalous) Ward identities 

D * JLk ~ ^T B A *J b B ulk - i eA R BC A - i £ aF A *J Mk = (179) 

D * Jtl - e [A A * Jfi = 
But once again, the fluxes into £ are non-trivial. These are easily computed 22 

^Ql = Cujj a (180) 

-\ab Ch f a a J> 



AQf = / e a Ae b (181) 



We now write the Ward identities in the edge for the covariant currents J^ ov = J^ ons +j a and 



' cov 



d*J cov = a H F (182) 
D * J a cov - i e «T b A * J b cov - i e aF A * J cov = („r (183) 

e [a A*Ji = C -fe a Ae b (184) 

Notice that the right-hand sides agree with the charge (in this case energy-momentum) 
entering the edge from the bulk. 

Let us now extend the above analysis to include curvature. One finds that the covariant 
anomalies, when written in terms of torsion and Levi-Civita curvature, become 

d*J CO v = a H F (185) 
D*J^ ov -i ea T b A*J b cm] -i e aFA*J cov = ( H T a + K H (e a AdR-RT a> ) (186) 

e [a A*J± = l - (Ch ~ khR) e a A e b (187) 

In the next section, we will show that it is possible to derive these identities from the edge 
point of view using the Fujikawa method with a suitable choice of regularization. 

Given the above results, it is natural to ask if the torsion terms in the diffeomorphism Ward 
identity could be removed by the addition of local counterterms. Indeed, a shift of the frame 
current 

* J a cov *J a cov -\{Ch~ kh R) e a (188) 



22 Here, we disregard the normal component AQ^- Since this is related to a bulk diffeomorphism normal 
to the edge, we expect that it is related to extrinsic rather than intrinsic edge properties. 
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would make it symmetric. This shift however does not come from a local counterterm in 
the boundary theory, as one can easily see after close inspection, nor is it an ordinary 
improvement term. 23 In fact, it amounts to shifting the bulk effective action by 

AS oddMlk = - \ I (h e A A T A + I [ k h Re A A T A (189) 

1 J M 1 J M 

One of our main precepts is that divergences that appear in the bulk are common to all 
phases. Thus shifting the values of (#, etc. by finite counterterms simultaneously in all 
phases is allowed, but this does not change the differences in their values between phases. 
Therefore, we cannot avoid having a torsional response in one of the two phases. It is this 
invariant information that is encoded in the covariant anomalies of the edge theory, and 
these are the important physical effects. 

It is also interesting to see from equations (186) and (187) (also see Eq 141) that the viscosity 
is shifted by a term proportional to the Ricci scalar so that if our space-time manifold is 
trivial in the time direction and has constant spatial curvature then the viscosity is shifted 
by an amount proportional to 1/R 2 where R is the radius of curvature. 



8 1+ld Anomalies from Fujikawa method 

We will now derive the covariant Ward identities discussed in the previous section, from the 
edge point of view. We will use standard methods that produce the covariant anomalies, 
and the novelty of the calculation is that we will produce the torsional contributions to the 
anomalies. In so doing, we will encounter divergences associated with the torsional terms. 
Our context provides these divergences with a clear interpretation, as the ultraviolet cutoff 
of the edge theory is determined by the mass gap in the bulk, and their presence is linked 
with bulk transport properties. 24 

The chiral fermions localized on a 1+1-d space-time manifold £ coming from the boundary 
of the manifold M couple to the frame e a on E. For simplicity, we will assume that the 
geometry near £ is separable, with a frame of the form e A = (Ndx,e a ). For our purposes, 
it will also suffice to ignore extrinsic couplings to the chiral fermions because these do not 

23 An improvement of the frame current is a current redefinition which makes it symmetric, but does not 
modify it's conservation equation. 

24 Similarly, in 3+1 dimensions, [34] torsional contributions to the chiral anomaly have been found, with 
divergent coefficients, although the interpretation of these divergences was not understood. These divergences 
can be given a physical interpretation in terms of 4+ld topological insulators. [35] 
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affect the covariant anomaly computations which are of interest here. 25 

Let us quickly review the Fujikawa method for computing covariant anomalies. Our dis- 
cussion here mainly follows [36, 31]. The main point of the Fujikawa method is that the 
variation of the fermion measure under symmetry transformations leads to anomalous Ward 
identities. For a Dirac fermion one defines the measure as follows: expand \I> and \I/ in 
terms of eigenfunctions $ m of a self adjoint operator, conventionally chosen to be the Dirac 
operator 

t^m = A m $ m (191) 

^ = ^a m $ m , * = J2b n ®n (192) 

m n 

and define the measure as [cM/d 1 ]/] = Yimn^mda n - However for a left-chiral fermion x/j, 
the operator T() L = $>|(1 — 7 5 ) is not self-adjoint. Thus, ip must be expanded in terms of 
eigenfunctions m of 'PlPl an d ^ must be expanded in terms of eigenfunctions Xn of 'Pl'Pl- 
Under a symmetry transformation T : ip — > ip' = ip + St"^, the measure could transform in 
general 

[dtfj'chp) = e-^ AT [d^cUp] (193) 

When this happens, the classical Ward identity gets modified by the anomalous correction 
At- Let us now study this in detail for diffeomorphisms, local Lorentz transformations, and 
Nieh-Yan-Weyl transformations. 



8.1 Diffeomorphisms 

Recall from (48) that under a covariant diffeomorphism generated by a vector field £, we 
have 

Sifj = VgV (194) 

For simplicity, we restrict ourselves to volume preserving diffeomorphisms. Then to linear 
order in £, it is easy to check that the measure transforms as 



'#'] = exp I - ^ / V 0h4>mC^i<Pm + ^2 V0 ^XnC^iXn ) 
\ m J n J 



[##] (195) 



2.-> 



The only extrinsic coupling to the chiral states is through a term in the effective action of the form 

S ext ~ ^- / ln(N)R volx (190) 



where R is the Ricci scalar on S. This term ensures that while the edge theory is anomalous under a Nieh- 
Yan-Weyl transformation of the frame on S, there is no anomaly due to a Nieh-Yan-Weyl transformation of 
the bulk frame. 
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Thus the corresponding Ward identity (50) gets modified to 

i [ ?(D* J a -i ea T b A*J b -i ea FA*J) = J2 [ vol s (0 m rV a m - XnCVaXn) = -Tr 7 5 rV a 

(196) 

Clearly, the trace is ill-defined by itself, and needs to be regulated. Customarily, it is 
regulated using the heat-kernel regularization in Euclidean space 

- Tr 7 5 fV/ 2/A2 (197) 

where A the ultraviolet cutoff is taken off to infinity. However for the edge theory we 
consider, the ultraviolet cutoff A is of order m, the bulk mass gap, since the spectrum of 
localized edge modes of the bulk Dirac operator only exists for energies E < |m|. This issue 
is irrelevant in the torsionless case, because the leading terms in the anomaly are finite and 
cutoff independent. However, in the presence of torsion we find a quadratic divergence in the 
anomaly if regulated naively. Moreover, the divergent term cannot be removed by a local 
counterterm. 

The guiding principle in choosing the appropriate regularization must then be bulk-boundary 
matching — the non-conservation of charge as manifested in the covariant anomaly must 
match the influx of charge due to the parity violating terms in the bulk action. Fortunately, 
a minor generalization of the results in [31] readily implies that (in a separable geometry) 
the flux of the bulk frame current in a given phase into the edge is given by 

- 1 1 
-g- WMi )^fV. (i ^ (198) 

where are the masses of the bulk fermions, including Pauli-Villars regulators. Notice that 
this is exactly the trace that was obtained in the Fujikawa formalism, albeit in a regulated 
form. Therefore, it is clear that we must regulate the Ward identity (196) as 

r - i i 

i / ?(D * J a - i ea T b A * J b - i ea F A *J) = ~AJ2 ^ sign(M,)Tr 7 5 £ a V a - 

(199) 

with Mq = m the mass gap in the bulk, and Mj for i = 1, • • • N being the masses of the 
Pauli-Villars regulator fermions in the bulk. The symbol A indicates that the anomaly is the 
difference between the flux from the non-trivial phase and the flux from the trivial phase. 



In order to compute the trace (in Euclidean space), we can rewrite it as 



N poo 

~ A J2 or(T\ CiMi / dt r 1/2Tr ^ V « e_i 



(-t> 2 +M?) 



(200) 
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The asymptotic expansion corresponding to this trace in 2 dimensions is given by (see Ap- 
pendix A) 

TV rt-V/" =, -i fe(±T.+ -± j +.Adk-± ; i i T.+ ...) (201) 

where the ellipsis denote terms higher order in t (which are unimportant here as they will 
give higher order terms suppressed by inverse powers of the cutoff). The integral over t 
in (200) diverges as e — > for the first term above, but the divergence is cancelled by the 
condition Y2iLo CiMi = on the regulator masses. Using the expressions for regulated bulk 
coefficients (148, 149, 150), we get the Ward identity 

(D * J a - i ea T b A * J b - i ea F A *J) = ( H T a + k h (e« A dR - R T a ) (202) 

where (h and Kh are the regulated coefficients in the non-trivial phase. Note that this is 
exactly what we found in our analysis in the previous section (see (186)). 



8.2 Local Lorentz transformations 



The change in the measure corresponding to Lorentz transformations Sip = |0 O 67 if) is given 
by 



'#'] = exp (^-^Tr 6 abl 5 r b) j (203) 

Following the discussion of the diffeomorphism anomaly in the previous section, we regulate 
the Ward identity as 



r 1^1 i 

^ el» A * J* = --A g - Ci sign(M ( )Ti Srfi* _ 



(204) 



Using the asymptotic expansion in 2 dimensions 



Tr 7 VV^ ~ / -%t ah ( - ^-R + ■■■ ) vols (205) 



E 



we find (187) 

e [a A * Jl = 1(Ch- k h R) e a A e b (206) 



43 



8.3 Nieh-Yan-Weyl transformations 



Under a Nieh-Yan-Weyl transformation Sip = — \uip-, the measure transforms as 26 

[dip' (hp) = exp QlY a J (207) 

The situation with Nieh-Yan-Weyl transformations is slightly different - the trace here cannot 
be interpreted in terms of the flux from bulk parity violating terms. However, there is no 
problem with regulating the anomalous Ward identity by the traditional heat kernel method 

i J a e a A * J a = ^Tr oe^'^ (208) 

This is because the quadratically divergent terms in this case can be removed by an appro- 
priate local counterterm. We will therefore only be interested in the finite and universal 
piece in the trace. Using the asymptotic expansion 



(. 



Tre '^Mi-2i^ + "') 1,ofe (209) 

A 2 

and shifting the current by a local counterterm J a — > J a — -^e a , we get 

e a A*J a = -— Rvoh (210) 

4o7T 

Indeed, this agrees with the well-known result T M M = — Z R, because for a left-handed 
Weyl fermion (cl,cr) = (1,0). 



9 Spectral Flow 

Given the close analogy between the Hall conductivity and Hall viscosity that we have 
discussed above, it seems reasonable that there ought to be physical adiabatic processes 
that induce spectral flow in the interface in each case. The case of the Hall conductivity is 
well-known, and we briefly review it here. 

26 It is important to note that both chiralities contribute with the same sign to the Nieh-Yan-Weyl anomaly, 
which is therefore proportional to (til +^_r)- On the other hand, the anomalies corresponding to diffeomor- 
phisms and local Lorentz transformations come with (til — nji). 
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9.1 Spectral Flow and Hall Conductivity 



We consider a gauge field on a spatial cylinder of length L in the x-direction and radius R 
in the ^-direction 

A = E y tdy (211) 

where E y is a constant. This is equivalent to 

F = EydtAdy (212) 

* 3 F = -E y dx (213) 

Thus we have a constant electric field in the y-direction which we imagine resulting from 
the threading of electromagnetic flux along the cylinder. We can parameterize E y = — 2w ^ RT 
where q is the charge and T is the time it takes to thread one flux quantum into the hole of 
the cylinder. 

Given the effective action for the bulk charge response S^/LA] = J M ^cr H A A dA), where 

2 

a H = — 2- and A is the electro- magnetic gauge field, the expectation value for the current is 
given by J = ou *3 F. Thus we find the bulk current response to the electric field is 

J = q \ H dx=—^—dx (214) 
2ithqRT 2ttRT v ' 

This means there is a constant current density in the ^-direction, and over the time T we 

build up charge 

AQ= dt dyJ x = q. (215) 
Jo Jo 

From the point of view of the intrinsic boundary theory, which consists of chiral fermions 
of one chirality or the other, this increase in charge is an anomalous process. In 1 + 1, the 
chiral symmetry is anomalous d* 2 J$ = 2^-^(2), where Fm) is the gauge curvature in 1 + 1. 
In the present case, this is the pull-back of the bulk field, which is just F( 2 ) = E y dt A dy. 
So the axial charge changes in this process by AQ 5 = L d * 2 J5 — — 1. This change occurs 
as a chiral fermion is pumped from one edge of the cylinder to the other through the bulk 
Chern-Simons response. 

We can get a simpler pictorial understanding of the anomaly by considering the energy 
spectrum of the chiral fermions. The Hamiltonians for the left- and right-handed chiral 
fermions are 

H R = v(p - qA) H L = -v(p - qA) (216) 
where the vector potential is A = Substituting this form into Eq. (216) we find 

fi>v ( t\ hv ( t\ . . 
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Figure 4: (a) Energy spectrum from Eq. (216) at time t = 0. Right/left handed fermion spec- 
tra are represented by positively/negatively sloped lines. The filled/empty circles represent 
occupied/unoccupied states, (b) Energy spectrum at time t = T where one flux quantum 
has been threaded through the spatial ring. The spectrum returns to itself but the state 
occupation changes. One electron has been added to the right movers, and one has been 
removed from the left movers. 



where n is an integer labeling the discrete momentum modes p = Assuming that T is 
very large so that the spectrum changes adiabatically, we find that the spectrum flows as time 
increases. At a time t = T, or in fact at any multiple of T, the spectrum returns to its initial 
configuration, yet the system as a whole has changed because the state occupation changes. 
When t = rT for integer r there have been r flux-quanta threaded into the circle on which 
the chiral fermions live. For each flux quanta threaded an electron is transferred from the 
left movers to the right movers as illustrated in Fig. 4. Thus we reproduce our calculations 
from above by observing the transfer of electrons during the spectral flow process. 



9.2 Spectral Flow and Hall Viscosity 

Now we would like to understand the momentum transport due to the Hall viscosity using a 
spectral- flow type argument similar to the case of charge transport. On the spatial cylinder 
of length L in the x-direction and radius R in the y-direction, consider the co-frame 

e° = dt, e 1 = dx, e 2 = (1 + h{t))dy (218) 
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Figure 5: (a)Energy spectra for left (blue) and right (red) handed chiral fermions from Eq. 
222 at t — (b) Energy spectra for t > assuming bi = &r = b < which gives an increase 
to the velocity for both branches of chiral fermions. Note that when compared to the t — 
case there are states that were occupied chiral fermion states that have been pushed past 
the cut-off scale and the states at p = are unchanged. During this process no states cross 
E = and there is not a conventional notion of spectral flow at low-energy. Both figures 
have the same momentum discretization spacing, but different velocities which leads to a 
different number of states within the cut-off window. 

where we will parameterize h(t) = 2 ^ T where T is a very large time-scale so that the change 
is adiabatic, and b has units of length. For simplicity, we will choose the connection uab = 0, 
for which the given frame is torsional. This configuration represents the threading of torsion 
flux T 2 = 2^rfdt Ady, i.e. a dislocation into the ring with a time-dependent Burgers' vector 
tangent to the ring with length bt/T at time t. To calculate the bulk energy-momentum 
flow we must introduce a covariant Killing vector field 27 £ = £ a e a = d y . From our previous 
discussion of Hall viscosity response, the energy-momentum flux along £ through a constant 
x slice T, x is given by 

L ( ' tJ ^ L ^ = i "L ( 1 + 2^ir) 2^Rf dth * (219) 

This leads to a transfer of momentum from one edge to the other through the bulk of 
the cylinder. From the point of view of the edge chiral fermions localized at x = and 
x = L, this is an anomalous process. For instance at x — 0, define the chiral momentum 

27 We call a vector field £ covariantly Killing if the co-frame is preserved under a covariant diffeomorphism 
along f , i.e. if D£ a + i s T a = 0. 



47 



Pi{t) = f t £ a * J a , where 7^ is the spatial circle at time t and J a the intrinsic 1+ld frame 
current on £ x =o. Then the anomalous conservation law (183) becomes (in the absence of 
U(l) gauge fields) 

dP L , f 2nR ( bt \ b J t bt \b 



dt 



f vti ( bt \ b , / bt \ b , . 

W. u 5if w*^ 1 + « f (220) 



Let us now understand the anomalous momentum transfer in terms of the spectra of the left 
and right handed chiral fermions localized at x = and x = L respectively, with co-frame 
fields parameterized by b^ and bn with the Hamiltonians 

hvk 
hvk 

h l = ~Y^r ( 222 ) 

where u is the chiral fermion velocity and we have assumed the Hamiltonian is acting on 
translationally invariant plane-wave states. We show the energy spectra at two different 
times in Fig. 5(a) and (b). In the figure we have indicated a high-energy cut-off governed 
by the scale |m|. This scale represents the energy at which the edge states of a topological 
insulator merge with bulk states and are no longer localized on the edge. In fact, for such 
topological insulators like the lattice Dirac model, the cut-off is exactly the insulating mass 
scale \m\. We have assumed that the energy states are filled up to energy E = as indicated 
by the filled circles in Fig. 5. The range of momenta that is occupied by right (left) movers 
is between p G ^f) , 0] (p € [0, f (l + ) • The total momenta of the right 

and movers at time t is 



v \ ~ IttRT J 

2*n . r 



m ( bf/t 
1 + 



v \ 2tcRT 



P t 0t) = ~2^KL ' Pdp h 



m ( bit 
1 + 



v V 2tcRT 



(223) 



As we have seen, the Hall viscosity is related to a stress-energy response and thus to the rate 
of change of momentum. We find 



P ( J ot) = -{ — ) [I + 



itl \ 2 / bit \ hb 



R " \hvJ V 2nRT J 2nT 
p(«) = (™\ a ( 1+ te\to* (224 ) 



\hvJ V 27rRT J 2nT 

We see that if we choose bi 7^ b^ then momentum is not conserved at all if we only consider 
the edge states and take into account transfers between the edges. Momentum of course is 
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still conserved globally because the excess/deficient amount of momentum gets trapped on 
some extra torsional flux that will appear in the gapped bulk region away from the edges 
when hi ^ 6r. For now we will fix b^ = bn = b to avoid this extra complication. 

Comparing equations (220) and (224), we see that the bulk and boundary momentum trans- 
port only matches for 



which is the same result we calculated earlier for the regulated Hall viscosity albeit with all 
the factors of h and velocity (speed of light) added back in. 

While at first glance it appears strange that the viscosity depends on the mass, we can 
clearly see the reason why this dependence is necessary by examining Fig. 5. The effect of 
threading a torsional flux (i.e. threading a dislocation) into the loop on which the chiral 
fermions propagate can be interpreted in one of two ways. Our choice of torsional flux (i.e. 
our specific choice of frame) means that if we travel around the loop at time t = T we enclose 
a Burgers' vector that is tangent to the ring of length b. Depending on the sign of b this 
implies that the ring looks either shorter or longer than its original length at t — 0. From this 
perspective we would think of chiral fermions with a fixed velocity but propagating on a ring 
with a time-dependent length (which will re-discretize the momentum modes as a function of 
time). The other interpretation is that the length stays fixed at 2nR but the chiral fermions 
are either traveling faster or slower depending on the sign of b. This is the interpretation 
represented in Fig. 5 where the velocity of the chiral fermions has increased at a later time 
but the momenta retain the original quantization scale. Thus we see that coupling to the 
U(l) electromagnetic field causes a translation in the spectrum, but the coupling to torsion 
causes a scaling of the spectrum. As a function of time the two chiral branches rotate in 
opposite directions around the fixed point where p = 0. This is because p = does not feel 
any effects of torsion since it is uncharged as far as torsion is concerned. So the torsional 
response is given by spectral scaling/rotation instead of spectral flow/translation. In terms 
of the discussion we used in the introduction this occurs because each momentum mode 
carries a different charge under torsion, while they all carry the same U(l) gauge charge. In 
fact, the state at p = does not even see the torsional flux and is unmodified since it carries 
zero torsional charge. 

9.3 Comments on Spectral Flow 

In this section we would like to discuss some comparisons between the anomalous electromag- 
netic and torsional responses in terms of the behavior of the edge state spectra. Anomalies 




(225) 
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Figure 6: Schematic illustration of the energy spectra for the case of (a) unmodified chiral 
fermions (b) insertion of an electromagnetic flux quantum (c) insertion of torsional flux. The 
red/blue circles represent occupied Right/Left handed chiral fermion states. The gray-dashed 
circles represent occupied states which are not chiral since they are beyond the energy cutoff 
at E = ±\m\ (represented by the dashed horizontal lines). Empty lines represent unoccupied 
modes. Note that in (b) none of the states are fixed and they all shift by the same amount 
whereas in (c) the states at p = are fixed and the states shift by amounts proportional to 
their momenta. 



connect low-energy physics with high-energy physics, and thus understanding both regimes 
is important for discerning the physics. For chiral fermions arising as the boundary theory of 
a higher dimensional massive Dirac model, both regimes can be understood. The low-energy 
physics is that of linearly dispersing chiral fermions that propagate along the 1+1-d bound- 
ary; this is the theory we have been carefully studying. In the high-energy regime there is 
an energy cutoff when E = ±|m| at which the chiral branch of states, which was localized 
on the boundary of the system, spreads into the higher-dimensional bulk and couples with 
the boundary states on the opposite edge. At this energy, and beyond, the fermions are 
no longer chiral, but in fact form Dirac fermions through their coupling with the opposite 
chirality edge states on the other boundary. Thus, states which lie outside of E = ±|m| are 
not chiral, and not localized on an edge. 

In Fig. 6 we compare the behavior of the spectra for both anomalous processes. Fig. 6a shows 
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the unmodified spectrum at t = for left and right handed branches while Fig. 6b, c show the 
modified spectra at a later time for the electromagnetic and torsional cases respectively. The 
red/blue colored circles represent occupied states, while the gray (dashed) circles represent 
occupied states beyond the cut-off which are not chiral. Empty lines represent unoccupied 
states. The modified spectrum for the electromagnetic case is simply a shift compared with 
the original spectrum, and is easy to understand since every state, no what the energy or 
momentum, shifts the exact same amount. The reason the shift is uniform over the entire 
spectrum is that all of the states carry the same charge q under the electromagnetic gauge 
field, and thus all couple minimally with the identical coupling constant. During the process 
of threading one flux quantum we see that one unoccupied state from the R branch passes 
beyond the cut-off and becomes non-chiral, while a formerly non-chiral occupied state in 
the R branch enters the edge state range and attains a right-handed chirality. That is to 
say there is an additional chiral fermion in the R branch after the process. The opposite 
phenomenon occurs for the L branch. Since the momenta of all states shift the same way, 
the notion of whether an anomalous charge appears or does not appear is independent of the 
width of the cut-off (though it does matter where the cutoff is centered around, which is to 
say the value of the chemical potential compared to the cut-off can be important). Thus we 
can easily imagining sending \m\ — > oo and we would still draw the same conclusion that one 
charge is produced for each inserted flux quantum. This even makes some sense (although 
it is imprecise) if we start with the cut-off at infinity since we know the states at the highest 
energies are shifted by the same amount as the states at the lowest energies since they all 
have the same electromagnetic charge. 

The torsional case is more complicated to understand. In this case the spectrum is stretched 
away from the states at p — 0. Unlike the electromagnetic case, where none of the states 
remain fixed when flux is inserted, under torsional flux insertion the states at p = do not 
change. Since each state is charged differently under torsion, i.e. the momenta of all the 
single-particle states are different, each state translates by a different amount. Of course, 
a momentum dependent translation that depends linearly on p is nothing but a scaling of 
the momentum, and thus a scaling of the energy since -Er/l — ±vp. We immediately run 
into a problem if we do not have a cut-off because the states at higher energies have higher 
charges and thus it is not clear how to interpret spectral changes if we allow p — > oo. For any 
well-defined system, such as a real material, or a properly regularized theory, this is not a 
problem because the momentum range over which chiral fermions exist will always be finite. 

As shown in Fig. 6c, when the torsional flux is inserted some states are scaled beyond the 
cut-off and are no longer localized, chiral modes. If we pushed the cut-off further out in 
energy there will be more and more states that get scaled beyond the cut-off per unit time 



51 



since states in a larger range of energy near ±\m\ will have momenta that are large enough 
to scale them beyond the cut-off. We see that while the same states stay occupied, and 
interestingly, no states cross E — 0, the amount of states which are actually localized on 
the edge is reduced. Since none of the states flow through E = 0, a connection to an index 
theorem, if one exists, must arise from a new mechanism. Note that if we considered the 
opposite sign of the Burgers' vector the states would be scaled downward in energy so that 
the number of states within the cutoff would be denser. Thus, the momentum transport 
does not come from a change in occupation of the states due to a spectrum shift, but instead 
from a decrease/increase in the density of chiral modes that lie within the cut-off. The 
electromagnetic modification to the spectrum thus acts like a rigid, incompressible flow of 
levels while the torsional modification allows for changes in the density of levels around some 
fixed point; in both cases the cut-off energy acts as a source/sink of levels. The change in the 
momentum of a chiral branch depends crucially on how many levels pass through the cutoff 
and what momenta those levels carry. Since both the number of states, and the amount of 
momenta the states carry, is linearly proportional to the cut-off energy \m\ we expect that 
the viscosity, which represents the rate of momentum change, would be proportional to \m\ 2 
as we have already found. 



10 Conclusions and Remarks 



In this paper, we have presented detailed calculations of physical quantities in Chern in- 
sulators, particularly in 2+1 dimensions. In this model, the effective action (generating 
functional) takes a Lorentz-invariant form. The leading parity-violating effect, the Hall vis- 
cosity, is encoded in this effective action through a term involving torsion. The interpretation 
of this term has been hampered by the appearance of power divergences. This should be 
compared with the (integer) quantum Hall systems, whereby similar computations give finite 
answers directly. This occurs simply because there is an effective cutoff set by the magnetic 
length, and the computations do not encounter divergences because there are a finite number 
of states (on the spatial torus, say). Nevertheless, careful management of the divergences 
that appear in the Chern insulator lead to a consistent picture in which the difference in 
Hall viscosities of distinct topological phases is independent of the cutoff and set by the mass 
gap. As one might well expect, this difference is also encoded in the structure of anomalies 
on an edge between such topological phases. Indeed, we have shown that the covariant dif- 
feomorphism anomaly of a 1+1 chiral fermion includes a torsional term that fully accounts 
for the jump in Hall viscosity across the edge, just as the chiral anomaly accounts for the 
jump in Hall conductivity. 
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The appearance of a chiral gravity action with negative cosmological constant in the lead- 
ing terms of the effective action is intriguing. Certainly this is consistent with the chiral 
nature of the edge theory, encoding cl — cr = 1. Is there some relation here with dynami- 
cal/holographic gravitational systems? In our context of encoding transport coefficients, the 
answer is surely negative - there is no role played by gravitational field equations or any 
particular solution thereof. However, given that dislocation/disclinations are analogous to 
vortices, effective dynamics may emerge as the result of a condensation mechanism [37] . We 
will report on this idea elsewhere [3 5]. 

One of the features that we have encountered above is the presence of an anti-symmetric 
part for the stress-energy tensor. We reiterate here that its presence is not forbidden by 
principle but is usually not considered simply by assumption of the nature of the matter 
being discussed. However, in any medium in which the microscopic degrees of freedom have 
local spin degrees of freedom, we can expect that such effects may be manifested. This would 
apply both to materials as well as presumably to hydrodynamic fluids. 

We expect that much of the physics that we have discussed in this paper will appear anal- 
ogously in higher dimensional systems as well, although the details will be quite different. 
We will return to such cases elsewhere [35]. 
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A Appendix: Asymptotic expansions 

On a number of occasions, we have encountered traces over Dirac fermions of the form 

Tr 7 V^, Tr 7 5 e a V a e^ 2 , Tr 7 5 7 a V^, Tre^ 2 (A.l) 

and in particular, their asymptotic expansions (in powers of t) in the limit t — > 0. We can 
use Af = 1 supersymmetric quantum mechanics to evaluate these expressions. 28 We will not 

28 The calculation of asymptotic expansions in presence of torsion has appeared in [38, 28, 39]. The use of 
J\f = 1 Supersymmetric quantum mechanics (SQM) in computing Chiral anomalies or Atiyah Singer index 
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provide details, but rather only sketch the essential ideas involved; see [36, 45, 46, 47] for 
details. 



Let E be a manifold with metric g^, and a torsional connection Wj. a & = £j^ ab + C^ ab . The 
action for M = 1 SQM in the presence of torsion is given by 

Ssqm = J dt (^Qijx^ + l - X a (S abX b + x k u k , abX b ) - i q j;X k x a X b H kab - ~N abcdX a X b X c X d 
+ %c[c + %x k A k c) + l -cF abX a X b ^j (A.2) 

where x % are local coordinates on E, X a are one-component real fermions, while c and c 
are one-component complex fermions. The theory is invariant under the supersymmetry 
transformations 5x l = ie X \ 5 X l = —ex 1 , with the supercharge 

Q = iX a eL(Pi ~ l^bcXV + cAc) - q j-±H a;bcX a X b X c (A.3) 

(pi being the momentum conjugate to x l ), and the Hamiltonian % = —Q 2 . Upon quantiza- 
tion, we must replace pi —> —idi and X a —> ^pfl a ■ The supercharge becomes Q = -^t> + •••■, 
while the Hamiltonian is H = —\Tp 2 + • • •, upto operator ordering ambiguities indicated by 
• ■ -. Further, the fermion number operator in SQM, (— is proportional to the chirality 
matrix 7 5 . This is what allows us to compute traces of the type (A.l) - the Hilbert space of 
M = 1 SQM essentially furnishes a representation of Dirac fermions on E. 

For instance, consider Tr 7 5 e*^ . In SQM, this is proportional to the Witten index Tr (— l) F e~^ 
with t — Such a trace over the Hilbert space is easiest to compute using the path inte- 
gral representation. To handle the operator ordering ambiguities, we follow the time-slicing 
prescription for the path integral [45], at the expense of the counterterms 

t - 1 Jir* r l _l 1 , Mt) r ,(iT) i ' ab ^ 1 u n a ; bc (k A\ 

Let - «1 jk + Vfc^ - ~16~3! ; ^ ' 

The path integral corresponding to Tr (— l) F e~ l3n is then given by 

Tr(-l) F e^= I [dx i d X a da i db i dc i )e-^ dtLB (A.5) 

J PBC 

where are commuting ghosts, hi and q are anti- commuting ghosts, 29 and Le is the Eu- 



densities on torsional backgrounds has been discussed before in [40, 41], and in the special case of vanishing 
Nieh-Yan four form in [42, 43, 44]. 

29 The ghosts are introduced to exponentiate factors of det(e) which arise due to insertion of complete set 
of position eigenstates in the discretized path integral. 
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clidean time Lagrangian given by 

Le = l 2 !J, r r'-r' + 5 W* 6 + ^wgfcXV + yW^xV 

+ c(c + i fc A fc c) - |cF a6X Yc + ^-(oV + 6V) + L ct (A.6) 

Here and a 1 have periodic boundary conditions, x a have periodic boundary conditions 
because of the (— 1) F in the trace (which is what the subscript PBC indicates), and 6j, Cj, c 
and c all have anti-periodic boundary conditions. In the absence of (— \ a acquire anti- 
periodic boundary conditions (APBC). Finally, the (3 — > limit is just the weak coupling 
limit in SQM, where we can do perturbation theory. In this way, M = 1 SQM allows us to 
compute the asymptotic expansions in (A.l) using standard techniques of field theory. For 
instance, using the method described above, we find the asymptotic expansion for Tr 7 5 e*^ 
in d = 4 is given by 

Tr 7 V^ ~ / (^T a AT a -R ab Ae a Ae b +^FAF + -^—R ( --^ a b AR ( --^ b a 
+ ^- 2 dd\T a AT a -R ab Ae a Ae b )+0{t)) (A.7) 

The same procedure can be applied for computing the other asymptotic expansions in (A.l). 
We state some of the results relevant to the calculations of sections 5 and 6. In d = 2 we 
have 

Tr 7 V^ ~ f (^-F + 0{t)) (A.9) 



2tt 



T r7 y'e*< H j^(±. r7 L k + 0(t) y,o h (A.10) 
H-^fV.e-* 2 c J^(J-T° + ^e"AdR-± r f)T° + 0(t)) (A.ll) 
while in d = 3 we need the result 



B Appendix: Results and Conventions 

We work in a mostly plus metric r] ab = diag(— 1, 1, . . . , 1) and take the Clifford algebra to be 

{y\ 7 fe } = 2r^ (B.l) 
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represented by numerical unitary matrices. 7 is then anti-Hermitian and all others are 
Hermitian. Furthermore, we have ip = ^7° . We use notation like jab = ^[la,1b], and we 
have 7°7 a (7°) _1 = — (7°)^ and ( 7 ) _1 ( 7 ab )t 7 ° = — 7 ab , etc. We note the following useful 
commutators and anti-commutators, valid in any dimension and signature 

{ 7 a , 7 bc } = 2 7 a6c (B.2) 

[ 7 a ,7 bc ] = 2(r/V-7/V) (B.3) 

{ 1 a \ 1 cd } = 2-f abcd + 2(77 bc r] ad - V bd r] ac ) (B.4) 

[ 7 ab ; 7 cd ] = 2(r] bc 'j ad - rj bd ^ ac - rj ac ^ bd + rj ad ^ bc ) (B.5) 

{ 7 a , 7 fecrf } = 2(r/ a V d + r] ac j db + V ad l bc ) (B.6) 

[ 1 \ 1 hcd ] = 27 a6cd (B.7) 
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